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Abstract 

We introduce two approaches to analyses of Borel sum rules. In the first method, we analyze 
the sum rules in the limit Af^ — s- oo. Here, the sum rules become very sensitive to the chosen 
threshold sq. We fix the threshold by setting the daughter sum rule equal to the meson mass in 
this limit. The second method introduces a Borel mass dependent threshold so{M'^). We choose 
functions so(M^) such, that the corresponding sum rule is not dependent on anymore. The 
relevant hadronic parameter is extracted at the most stable function so{M'^). The two methods 
differ in the notion of quark-hadron duality. Whereas the first method emphasizes errors coming 
from the simple duality ansatz, the second method is constructed such as to extract hadronic 
properties from the region where the duality approximation is fulfilled best. We use these 
modifications of the sum rule approach to extract values for the B-meson decay constant Jb, 
the semi-leptonic form factor f~^{0) at zero momentum transfer and the strong coupling g^,g^. 
The results coincide with other sum rule analyses. 
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Chapter 1 

Introduction 



In 1979 M. Shifman, A. Vainshtein and V. Zakharov proposed the so called QCD sum rules 
^ |2j - a method to relate hadronic properties with QCD calculations, where non-perturbative 
effects are systematically included. With these sum rules, from now on referred to as SVZ sum 
rules, one is able to express static hadronic observables like masses and widths in terms of a 
few universal phenomenological parameters, the vacuum condensates of quark and gluon fields. 
Today, sum rules and their various derivatives are by no means the only approach to predict 
hadronic parameters. Other methods are calculations in heavy quark effective theories and chi- 
ral perturbation theory. Compared to sum rules they can only be applied to limited regions of 
hadronic physics - systems involving heavy quarks and only light quarks, respectively. Further 
methods are lattice calculations and various quark models. The advantages and disadvantages 
of lattice results are discussed below. 

The basic idea of the SVZ sum rules is to represent hadrons by currents of quark fields, that have 
the same quantum numbers as the corresponding hadrons, and to construct the two-point corre- 
lator of a current and its hermitian conjugate. This structure is then evaluated in the euclidean 
region for large spacelike momenta, where the high virtualities of the quarks ensure that the 
main contribution to the two-point function can be extracted from perturbative calculations. 
Corrections are added by use of Wilson's operator product expansion (OPE), separating the 
hard and soft regime of the correlation function to take large distance dynamics into account. 

2 

This gives rise to a truncated power series in where is the scale introduced in the OPE 
to separate short and long distance dynamics. By the use of dispersion relations the two-point 
function is then related to the hadronic spectral density. Finally, a so called Borel transforma- 
tion is applied to improve the convergence of the power series on one side and suppress higher 
states and continuum contribution of the hadronic spectral density on the other side. The SVZ 
sum rules will be derived and discussed more intensively in section ITT] of this thesis. 

Over the last 25 years, the initial SVZ sum rules have further been developed and improved. 
They have been applied to three-point functions, baryonic currents and external fields and were 
used to determine the masses of the light quarks, decay constants, form factors, couplings and 
magnetic moments of mesons and baryons. They were extended to investigate properties of 
hadronic matter at high temperatures and densities and many more systems (see |31 E] for some 
examples and applications of various sum rules). Besides the 2-point SVZ sum rules, we will 
make use of the so called light-cone sum rules (LCSR) in this thesis. Here, the relevant non-local 
matrix elements will be parameterized in form of wave functions instead of the expansion in local 
vacuum condensates. This method turned out to be quite useful to extract hadronic properties 
from reactions including three particles. We will derive the LCSR in section 031 
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One of the major advantage of the sum rule approach is its flexibihty. Sum rules can be applied 
to a variety of problems and give reliable predictions of hadronic parameters. Compared to 
lattice results, which became more and more important in the last years, sum rules do not give 
very accurate results - usual uncertainties lie in the range of 10% to 20%. Nowadays, lattice 
results have smaller errors and will certainly achieve higher accuracy in the near future, when 
unquenched calculations will be feasible within an acceptable period of time. But until now, 
lattice calculations still include relatively high masses for the light quarks, which are related to 
the length of the lattice. Thus, the results have to be extrapolated to the small masses of the 
light quarks to give physical predictions, which reduces the accuracy of lattice calculations. Fur- 
thermore, they only give brute force numerical results, where the fundamental physical dynamics 
are hidden beyond the extracted numbers and barely observable. On the other side, the simple 
analytical sum rule method allows physical insight into the systems they are applied to and an 
interpretation of the results - at the expense of accurate numerical results. Furthermore, the 
uncertainties of the results can be traced back to uncertainties of several ingredients of the sum 
rules, like the error coming from the truncation of the condensate expansion of the correlation 
function and the assumption of quark-hadron duality, which is necessary to compensate for the 
scarce knowledge of the hadronic spectral density. Another advantage of the sum rule approach 
is the fact, that it can be continued to Minkowski space in an analytical way, whereas this con- 
tinuation of the numerical lattice results obtained in euclidean space sometimes is problematic. 
Thus, although they suffer from some inaccuracies, sum rules are still an important method of 
calculating hadronic properties. 

In this thesis we will propose two modifications to the usual sum rule calculations of the decay 
constant /b, the form factor /^(O) at zero momentum transfer q'^ = 0, and the coupling q^'b-k 
the heavy-light B-meson. To have an precise measurement, or at least an accurate prediction of 
these parameters is of high importance. The weak decay constant Jb contributes to calculations 
of B — 5- mixing amplitudes and is needed in determinations of the CKM-matrix element Vub, 
when it is extracted from leptonic decays of the B-meson. It contributes as well to the sum 
rule calculations of the form factor f^{q^) and the coupling gB*B-K- '^^^ form factor can also 
be used in predictions of the CKM element. Knowledge of the fundamental parameter V^b is of 
greatest importance in today's physics, since it contributes to the question of unitarity of the 
CKM-matrix and thus tests the standard model. 

This thesis is structured as follows. In the next chapter, we will introduce the basic methods 
and theoretical tools. Different sum rules and their ingredients will be discussed. Furthermore, 
some of the crucial points of the sum rules will be stressed, like the dual relation of the hadronic 
spectral function to the perturbative QCD calculation. 

In chapter |31 we will start with briefly reviewing the sum rule analysis for the weak decay con- 
stant fs- We will basically follow an approach from Jamin and LangeiS , using the running 
b-quark mass in the MS-scheme. Compared to analyses using the pole mass of the quark, this 
shows an improved convergence of the perturbative expansion in the strong coupling Us- In 
most sum rule approaches, the hadronic parameter is extracted at rather ambiguous values of 
the squared Borel mass and the threshold sq. In the sections 13.21 and we will introduce 
two new methods to extract a distinct value of the hadronic parameter. 

In section we will evaluate the sum rule in the limit — > oo. This is related to an approach 
proposed several times by Radyushkin|2| El- Hereby, uncertainties coming from the truncation 
of the expansion in the condensates essentially vanish and the threshold parameter sq can be 
fixed to one distinct value, thus leading to a unique numerical result for fB- The drawback of 
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this method hes in the emphasis of uncertainties coming from the simple duahty ansatz. 

Section 13.31 introduces another possible modification. Here, we will introduce a Borel mass de- 
pendent threshold so{M'^) instead of a constant sq. With this at hand, one can find functions 
so{M'^), for which the sum rule is constant in a given window of the Borel mass. In this case, 
also the ratio of the sum rule and its first derivative, which is often used to get a reading point 
(M^, So), is constant and equal to the squared mass of the ground state of the hadronic side. For 
each value of the sum rule, we find a corresponding function sq^M"^). We will give arguments 
to extract the value of the hadronic parameter, for which the function so(M^) is most stable. 

In the chapters |1] and El we will apply the two methods to the analyses of the B-meson form 
factor /^(O) at zero momentum transfer = and the strong coupling g^,^^. The LCSR 
of these parameters are proportional to the weak decay constant Jb and, in the case of the 
coupling, also to the decay constant of the i?*-meson fs* ■ We expect some cancellations of the 
intrinsic errors of the two methods, when the ratio with the decay constants, analyzed within 
the same methods, is taken. Our extracted results are close to the values obtained in other sum 
rule approaches. Chapter El will conclude. 
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Chapter 2 

Theoretical Foundations 



Quantum Chromodynamics is described by the SU(3) Yang-Mills Lagrangian: 

C-QCD = -\{G^,uf + Y^i^fiip- mf)ijf, (2.1) 

/ 

where G^i, is in the adjoint representation of the gauge group and ipf in the fundamental - 
/ indicating the different quark- flavors. Although Cqcd is made out of quarks, rather than 
hadrons, it is believed to incorporate all hadronic physics. However, due to confinement, indi- 
vidual quarks are never observed experimentally, and thus perturbative calculations, which start 
from this Lagrangian are limited. Calculations involving the free quark propagator are clearly 
only valid far off-shell, which is one of the features of the SVZ sum rules, where the perturbative 
coefficients of the operator product expansion are evaluated at large virtualities of the quarks. 

To lower scales the running coupling as{^) increases and the perturbative expansion breaks 
down. It even reaches a pole at a scale of order Aqcd ~ 200 MeV and perturbative expressions 
are only valid down to energies of about 1 GeV. With increasing separation of two color charges, 
the interaction becomes so strong, that the potential between the charges cannot be described 
by the exchange of single gluons anymore. A whole cloud of gluons will develop between the 
charges and one has to describe physics at this scales with methods of (chromo-) electric and 
magnetic fields. A perturbative approach is not feasible in this region. 

However, for the treatment of hadronic parameters perturbative QCD is still applicable, since 
the average separations of the color-charged partons in the hadrons are not too large. One can 
treat them perturbatively and add corrections to the expressions by considering interactions of 
the partons with the QCD-vacuum. The QCD-vacuum is described by slowly fluctuating fields 
and can be included in the calculations by introduction of the so called condensates. With the 
sum rules, a systematic treatment was introduced by Shifman, Vainshtein and Zakharov in 1979. 
The SVZ sum rules will be derived in the next section. 



Throughout this thesis we will use the MS'-renormalization scheme with dimensional regular- 
ization in perturbative calculations. Adopted expressions from other authors are also obtained 
using this scheme. The coefficients of the corresponding Gell-Mann-Low function for the running 
coupling are listed in appendix Being a mass-independent scheme - the /3-function is not a 
function of the quark masses - the application of this scheme leads to difficulties in treating the 
heavy quarks, since they do not decouple in non-physical quantities. If one tries to express the 
running coupling constant at a certain scale as a function of the coupling at another scale, one 
has to apply matching conditions if the mass-threshold of a heavy quark lies between the two 
scales. This is also described in appendix IaI 
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2.1 The SVZ sum rules 



The basic object of the SVZ sum rules is the two-point correlation function of a gauge invariant 
quark current and its hermitian conjugate: 

I^AB{q) = i I d^xe'^^0\T{jA{x),j^^iO)}\0). (2.2) 

Here, jA{x) = ip f{x)TA'iljf'{x) is a quark field bilinear with Lorentz structure and flavor content 
chosen such, that it reflects the quantum numbers of the interpolated hadron. The currents 
are injected into the vacuum to avoid long distance dynamics, which would occur in case of an 
initial or final hadronic state. A typical process, which is relevant for investigation of properties 
of the p- meson, is shown in figure ITTl Here, the T-product of two vector currents is considered. 
In this case, due to current conservation dfj^j'^ = 0, one can explicitly factor out the kinematical 
transverse Lorentz structure: 

^(.uiq) = {q^^qu - gV)n(g')- (2.3) 

In general one can always decompose the initial n-point function into various terms, where the 
kinematical structure is factored out, and then proceed with the objects containing the dynamics. 

The idea of the sum rules is to evaluate the correlation function in two different ways. On one 
side, H(g^) is calculated at large negative momentum transfer = —q^ ^ Aqcd, where it 
can be shown to be dominated by short distance dynamics This will be done within the 
framework of Wilson's operator product expansion (OPE), which gives additional corrections to 
the pure perturbative part in form of an expansion in vacuum condensates. On the other side 
a complete set of hadronic states will be inserted in the correlator to obtain hadronic matrix 
elements. Usually, the sum rules are setup to extract values for these matrix elements or hadronic 
properties which they contain, like decay or coupling constants. The two sides are related to 
each other by the use of dispersion relations. A Borel transformation is applied to both sides 
to improve the stability. It suppresses higher states on the hadronic side and removes arbitrary 
subtraction terms that might appear in the dispersion relations. In the following we will derive 
the sum rules for the B-meson decay constant fs- 

2.1.1 The Hadronic Side of the Sum Rules 

Starting point is the two-point function of two pseudoscalar currents of a bottom and a massless 
quark: 

n(g2) = i j Ae^5"(0|r{g(x)f75Kx),K0)^75g(0)}|0). (2.4) 

It was shown by Kallen and Lehmann that any two-point correlator n(g^) is an analytic function 
in the complex g^-plane with possible singularities and a branch cut on the real positive g^-axis. 
The singularities and the cut appear at momentum transfers g^, where intermediate states go 
on-shell. For negative real < 0, ^{q^) is always real and fulfills Schwartz' reflection principle: 

n(z)=H*(z*). (2.5) 

This identity can be analytically continued to the whole complex g^-plane. With this at hand we 
can represent n((5^), where = —tf', by a dispersion integral, relating the two-point function 
at negative momentum transfer to its imaginary part at positive (f'. We consider the integral 
over the contour depicted in figure [2?2l 
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Figure 2.1: Appearance of the quark-loop in e e — > e e -scattering. At large negative values of the 
momentum of the virtual photon, the quarks can be represented by their free propagator. 



c 

If n(z) vanishes sufficiently fast at \z\ — > cxd, the integration over the circle does not contribute 
and we are left with 



oo 



2^ 1 /■ n(s ie) - n(s - ie) 







With use of ()2.5() the numerator can be replaced by the imaginary part of the correlation function 
and we get the dispersion relation: 



oo 



IT J s + Q^ J s + Q^ -le 







In most cases the imaginary part does not vanish in the limit s ^ oo and one has to introduce 
subtraction terms to make the integral finite. In our case the two-point function has a mass 
dimension of two and an asymptotic behavior of n(s) oc in the limit s oo. We therefore 
setup the dispersion relation for n(g^)/g^ instead of n(g^): 



oo 





Instead of (|2.8j) . we now have an additional polynomial in rendering the initial integral finite. 
The application of a Borel transformation (see section I2.1.3P will remove the subtraction terms 
and exponentially suppress higher contributions to the spectral function. 

We can get an expression for the imaginary part of the correlator n(g^) by use of the optical 
theorem. It relates the imaginary part of a forward scattering amplitude with the sum over all 
intermediate states - and thus with the total cross section for particle production of these states. 
Inserting in (|2.4j) a complete set of states yields 

2Imn(g2) = Y,\ dr„(27r)45(^)(g -p„)(0|gi7561n)(n|6i75g|0). (2.10) 

Here, |n) is a single-/multiparticle state with quantum numbers of the quark current and 
indicates the integration over the corresponding phase space. In our case, the B-meson is 
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n(Q2) 


9^ 











Figure 2.2: Contour of the integral relating Il{Q^) at negative momentum transfers to its imaginary 
part at positive g^. The circle is taken to infinity. The dots and the dashed line at the positive g^-axis 
indicate bound states and multi-particle states. 



the lowest lying state. The corresponding meson-to-vacuum matrix element defines the decay- 
constant fs, 

m^,b\B) = (2.11) 

rub 

where nib is the b-quark mass and niB the B-meson mass. Using this definition and integrating 
out the phase space, we get 

1 m'^ 

-lmU(q^) = -^4^5{q^ - ml) + higher states. (2.12) 
vr 

Here, higher states stands for the contribution of higher lying resonances (radial excitations) of 
the B-meson and the continuum contribution of multi-particle states. Very little is known of the 
spectral density, so for the time being, we write 



i/mn(g2) = ^5iq' - m|) + /(g^)e(4 - q% (2.13) 
vr 

where the step function is explicitly written out to emphasize that the onset of the contribution 
of the higher states and the continuum p^{q^) is at Sq > m^. This ansatz, first resonance plus 
continuum, for the hadronic spectral density is typically used in most applications of sum rules. 
Only the lowest lying state enters explicitly the calculations, whereas the hadronic continuum will 
be approximated by the perturbative results (see section [2^ . Plugging (|2.13j) in the dispersion 
relation 1)2. 8() . we get 

oo 

where the ellipses stand for subtraction terms necessary to make the integral finite. 

Equation 1)2. 14() is the hadronic side of the sum rules. On the other side, the correlation function 
n((5^) will be derived by means of perturbative QCD. In the next section we will calculate 
the simple quark loop and introduce corrections to it, coming from an expansion in vacuum 
condensates. 
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2.1.2 The QCD-Side of the Sum Rules 

We start with calculating the contribution of the simple quark- loop (figure E^la) • This means, 
we contract the quark fields in the correlator and replace them by the free quark propagators. 
In the chiral limit of massless u, d-quarks, we get 

.n(,^) ^ {,,,^J^^,,,^l±0.} (...) 

times an overall delta function. The '3' arises from summation over color indices. Instead of 
calculating this amplitude directly, we again setup a dispersion relation. By inspecting the poles 
in (|2.15j) . it can be seen that the correlator n(g^) develops an imaginary part from > on, 
and we can write: 

oo 
ml 

We calculate the imaginary part of ()2.15() by applying the Cutcosky cutting rules to it. We get 
the discontinuity across the branch cut on the positive g^-axis by replacement of the propagators 
with delta- functions^ . 

Discll{q^) = -—J^{k' - qk){27ri)5{{q - kf){27ri)S{e - ml) (2.17) 
After integration and using Discll{q^) = 2ilmll{q'^) we get 



3 {ml - q^f 
8tt g2 



Imlliq') = iz ' (2-18) 



Inserting this expression in (|2.16|) the integral becomes divergent. Again, this will be cured later 
by application of the Borel transformation. 

First order corrections in the strong coupling to the quark- loop (figure IT^ b-d) are well known 
and we adopt the expressions from P|: 

A/mn(s) = ^M^ci(,), (2.19) 
27r^ s 



with 



c\s) = - + 2L?2 +log — log 2 + olor ^ 



b ~ ""b ^ ~ "% 

2 ^ _ ^2 _2 



+ log 2 + log ^ + log — (2.20) 

s — s s — 



Here, the dilogarithmic function Li2{x) satisfies 



U,(x) = - / ^^^^dt. (2.21) 



0(a3)-corrections to the heavy-light system were calculated recently by Chetyrkin and SteinhauserjOj 
in a semi-analytic way. These three-loop corrections are available as a mathematica package. 



^This is a purely mathematical treatment to calculate the imaginary part of 11 (g^) - putting the quarks on-shell 
is certainly not a physical procedure. 
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Figure 2.3: The free quark-loop (a) and 0(Q!s)-corrections to it (b-d). 



Up to now, the two-point correlator was calculated only within perturbation theory. One can- 
not expect to reproduce hadronic properties by only considering free quark propagators in the 
calculations. At some point perturbative QCD breaks down and non-perturbative objects have 
to be taken into account. This will be illustrated in the next section. 



Appearance of non-perturbative Contributions 

If one considers higher order corrections to the quark loop, the appearance of non-perturbative 
corrections can be seen. We illustrate this by considering corrections to the simple loop of 
massless quarks appearing in figure l2?H following Shifman|lf)j. The Adler function is defined as 
the first derivative of the two-point function: 

D{Q') = -^^'Q'^^- (2.22) 

With this definition the loop without a^-corrections is unity and the inclusion of a gluon- 
exchange gives: 

oo 

AD(Q2) = as{Q^) J dk^ Fik^,Q^) = ^. (2.23) 



In higher orders one has to replace the coupling by the running coupling as{k'^) in the 
integral. This replacement takes the increase of the strength of the coupling constant into 
account, if only small momenta k"^ are exchanged by the gluon. It is equivalent to replace the 
free gluon propagator by the full one, carrying the effective charge (to one- loop accuracy): 

a (k^) - (2 24) 

i + -^if- Oo iog 

The resulting feynman diagrams are often referred to as bubble-chain graphs. The insertion of 
the effective charge leads to a resummation of graphs of this type (see figure |23J- Inserting the 
running coupling, the integrand becomes infinite at the Landau pole /c^ = A^: 



47r 



A^ = Q^e "0"" . (2.25) 
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Figure 2.4: Replacement of the free gluon propagator by the effective bubble chain. 



Dividing the integral in two parts at a scale fj? ^ A'^, the IR-part can be separated: 

AD{Q')jn = a,(Q2) [ dk" F{k\ Q^). (2.26) 

J 1 + 6n loe ^ 



There is no prescription for how to skip the pole, thus the integral cannot be solved in an 

2 fc2 



unambiguous way. However, using the small /e^-behavior of F{k'^,Q'^) f ^! one can deduce 



the scaling-behavior 



AD{Q^)jR a ^. (2.27) 



Thus, the low momentum contribution to the Adler function scales like with an ambiguous 
numerical coefficient. However, since the Adler function is an observable, this term has to be 
cancelled by some non-perturbative object of mass-dimension four. In the operator product 
expansion, which will be introduced in the following, such non-perturbative effects are system- 
atically taken into account by an expansion in vacuum condensates. The vacuum expectation 
value of the squared field strength, (0|^G'^'^G^iy|0), may be seen as the corresponding object 
incorporating the non-perturbative effects of the bubble chain. 



Wilson Operator Product Expansion 

The idea of Wilson's operator product expansionjllj is to factor out the long-distance part 
below the scale ^ into the non-perturbative vacuum expectation values of gauge and Lorentz 
invariant local operators, the vacuum condensates. These are accompanied by /x-dependent 
coefficients, which can be calculated perturbatively. This means that any correlation function 
will be split into a perturbative expansion, taking into account the coulomb interactions and an 
expansion in condensates, which bear the long distance interactions of the quark and gluon fields. 

The effect of the product of two operators separated by a small distance x can be described as 
the effect of a local operator of the same quantum numbers. This local operator creates the 
same disturbance as the product and can be written in form of a linear combination in some 
basis: 

O{x)O{0)=Y,Cn{x)On{0). (2.28) 

n 

The resulting operators On are purely local ones, whereas the coefficients C„(x), which are c- 
numbers, now depend on the separation x. 
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In the case of the correlator of two quark currents, we have to expand the product of the currents 
into a series of local operators: 

r(x)rt(o) = ^Cr(x)0„(0). (2.29) 

n 

Since we are interested in the vacuum expectation value of the currents, only the local operators 
On which are Lorentz and gauge invariant, contribute to the expansion. Taking the Fourier 
transform we can write: 

i y"Ae'«-(0|r{i'^(x),rt(0)}|0) = Cr(g2).i^^A-(^2)^-^^^^M-^(^2)^^^^^^M-^ + ... . (2.30) 

Here, 1 stands for the unity operator, so the coefficient C^'^{q^) is just the perturbative result 
of the quark loop. The condensates do not depend on the momentum transfer and also do 
not depend on quantum numbers or any other structure of the correlator. They are universal 
non-perturbative objects. The two condensates which are written explicitly in the expansion are 
the quark condensate of mass dimension three and the four dimensional gluon condensate. The 
quark condensate is the order parameter of chiral symmetry breaking. It is always accompanied 
by a quark mass and therefore effectively a dimension four object. The ellipses stand for higher 
dimensional condensates. The dimensions of the coefficients can be extracted directly from di- 
mensional analysis. If the kinematical structure is factored out - in case of two conserved vector 
currents this is just the transverse Lorentz structure {qfj,qy — q'^g^u) - it becomes clear that the 
condensates are suppressed by {Q'^)~^^'^, where d is the dimension of the relevant condensate. 
Thus the expansion in the dimensions of the condensates is suitable, if the correlator is evaluated 
at large negative momentum transfer. 

The correlation function n(g^) does not depend on any renormalization scale n, so the expan- 
sion on the r.h.s has to obey the Callan-Symanzik equation. We introduce a normalization 
point fj, such, that it separates the hard and the soft regime of the correlator. The soft modes of 
the vacuum fluctuations are then described by the condensates, which become scale-dependent, 
whereas the hard modes are calculated perturbatively. One chooses a rather safe scale, where 
the coupling as{fi) is somewhat smaller than unity, so that the perturbative expansion and the 
calculation of the Wilson coefficients is valid. 

The coefficients of the condensates are calculated by cutting the lines of the low momentum 
quark or gluon. Thereby, the momentum of this particle is set to zero and the fields are treated 
as external. In the case we are interested in, the heavy-light B-meson, the relevant diagrams 
up to dimension 6 are depicted in figure 12.51 The crosses indicate the insertion of the vacuum 
condensates. 

In our case of a pseudoscalar current of one bottom and a massless quark we have to start with 
(|2.4j) . To calculate the Wilson coefficient of the quark condensate {qq), Cgq{q'^), we contract 
the b-quarks in the T-product and replace it with the free quark propagator, whereas the light 
quarks are treated as external fields (see figure [231 a) . We expand the quark-field q{x) around 
zero, to get an expression in the local condensates: 

qix)=qiO)+x''Df,q{0) + --- (2.31) 

q{x) = m + mD^^xf" + ■■■. (2.32) 

The ellipses stand for higher derivatives which contribute to condensates of higher dimensions. 
The second term vanishes by the equations of motion for the dirac fields and we are left with: 
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Figure 2.5: Relevant Feynman diagrams for calculation of the coefficients of the condensates up to 
dimension 6. Appearance of the quark condensate (a), the gluon condensate (b,c), the mixed quark- 
gluon condensate (d) and the four-quark condensate (e). 



The term proportional to {0\q{x)^q{0)\0) gives contributions of order 0{mq), when the expansion 
(|2..S2j) is inserted and we neglect its contribution. Evaluating the integral we get 

CUq') = (2-34) 

Besides the four-dimensional condensates in (|2.3U|) . there are further contributions to the sum 
rules from the condensates depicted in figure [T31 d and I2.5l e. The dimension five mixed quark- 
gluon condensate is usually parameterized by: 

{9sq(T''''YGU) = ^o{m)- (2.35) 

The quark-gluon condensate appears, if one expands the non-local quark condensate. Using 
Fock-Schwinger gauge, Xfj,A>^(x) = 0, to make the path ordered integral between the two quark 
fields unity, one gets an expansion of quark condensates with covariant derivatives between 
them. The condensate (qD'^q) is directly related to the mixed quark-gluon condensate. Thus 
rriQ stands for the ratio {qD'^q)/ {qq), which can be interpreted as average virtuality of the quark 
fields. 

The four quark condensates can be factorized into two two-quark condensates in the limit of 
large number of colors (corrections come with 7^): 

(04g) = asiqqf. (2.36) 

The calculation of the coefficients of the higher dimensional condensates for the heavy-light 
pseudoscalar currents can be found in In ^Bji Jamin and Miinz calculated the coefficients 
to the condensates up to dimension 5 to all orders in the quark masses. It is found that long 
distance contributions in form of mass logarithms of the light quark mass log m^//i^, which 
appear in the calculations starting at order n = 3, can be absorbed in the condensates, if non- 
normal-ordered condensates are used in the operator product expansion. This fact was later 
used to obtain a first order correction in the strong coupling to the quark condensate [S]. 

Throughout this thesis, we neglect light quark masses, and we adopt the results from Aliev |14j 
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for the coefficients contributing to the heavy-Hght correlator: 



1 1 



CGG{q^) = -^^r^ (2-37) 
12 — mj^ 



2n IGtt 1 ml mt 



^'"^"^ ^ " 27 (g2 - ^2)2 (^^ + 2(g2 _ ml) 2(g2 - ^2)2 j • (^-^^^ 

Gathering all the coefficients of the expansion we get for the QCD-part of the correlation function 
n(Q2) up to dimension 6 in the condensates: 



oo 



um= It — 7^^^ 



(S + Q2) (m2 + Q2) 



/- \ , '^b / 1 \ 2/- \ , 



IGtt m2 > 2 



27(m2 + Q2)\^ 2(m2 + Q2) 2(m2 + Q2)2 



(2.40) 



where 



and C"'^(s) is given in (|2.2()|) . 

The values of the condensates do not depend on the channel one is investigating. We take the 
following typical numerical values for the condensates. 

The value of the quark condensate can be extracted from the Gell-Mann-Oakes-Renner relation 
obtained by current algebra: 

(m„ + ma) {uu + dd) = -mlfl (2.42) 

and we take from |15j : 

{qq){2 GeV) = -(267 ± 16)^ MeV^. (2.43) 

The gluon condensate is extracted from sum rules themselves, thus one channel is sacrificed. 
Instead of giving a prediction to a hadronic matrix element, this value is taken by experiment 
and the condensate adjusted such, that the sum rules reproduce it. From charmonium sum rules 
we get: 

(0|— G'^'^G.^IO) = 0.012 GeV^ (2.44) 
vr 



This value and the next two parameters are taken from BelvaevjlBj: 

ml{lGeV) = (0.8±0.2) Gey2 (2.45) 
as{qq? = 8 • 10"^ Ge^^ (2.46) 

We did not quote an error on the gluon and the four quark condensate. We will later take an 
uncertainty of 100% on these values, which will still be numerically negligible. 

The operator product expansion is an expansion in the soft non-perturbative condensates. This 
power series is only valid up to a certain dimension, where contribution from direct instantons 
appear. These short distance vacuum fluctuations give rise to an additional exponential term 
and the expansion in the condensates breaks down. However, the instanton contribution is 
numerically negligible. 
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2.1.3 Equating the two Sides 



The Borel Transformation 

After evaluating the two-point correlator in two different ways, the phenomenological, hadronic 
part and the QCD-part are equaled. To improve the stability, Shifman, Vainshtein and Za- 
kharov proposed a so called Borel transformation. Applying it to the two parts will remove the 
subtraction terms of the dispersion relations and improve the convergence of the power-series 
expansion at the same time. The algebraic operator is that of an inverse Laplace transform: 



B = lim lim 



m^(_^. (2 47) 



Applying this operator on any negative power of {s + Q ), we get 

After the Borel transformation, for sufficiently low values of the squared Borel mass the 



integral over the imaginary part of the correlator (|2.16|) becomes finite and also the contribution 
of higher states in ()2.14|) gets suppressed. The truncated power series of the contribution of 
the condensates becomes an expansion in negative powers of times an exponential function, 
which breaks down for — > 0. This gives a lower bound on M^. On the other hand, the 
higher the Borel mass is chosen, the more the uncertainties in the hadronic spectral density are 
emphasized. Thus, in the typical sum rule applications a working window of the squared Borel 
mass is defined. The boundaries are taken such, that the contributions from higher states 
and from the condensate expansion, respectively, do not exceed a certain percentage value of 
the whole expression. One hopes that inside this window the corresponding hadronic parameter 
the sum rule is set up for, is close to reality. The issue of which value of the Borel parameter 
should be taken as reading point, will be discussed several times throughout this thesis. 



The final Sum Rule for the B-meson Decay Constant 

Borel tranforming (|2.14|) and the expansion in local operators (I2.4U|) . we get after equating both 
sides: 



m%,fl, _!!Lb 
mt 



CO 

j dsp 



+ 



dsp^{s)e + 



, , I ,as m, 

-mbiqq) + —( — GG) ^ 

12 n ' 2M2 



IGtt 1 
' 27 M2 



1 



mt 



ml. 



4M2 12M4 



1_ "^b 



a. 



2M2 J 

2 



ml{qq)- 



e 



(2.49) 



where p^^^{s) is given in (|2.41j) . The hadronic spectral function p^{s) is not known. By the 
assumption of quark-hadron duality, it is related to the perturbative part of the correlator. It 
is assumed that at large momentum transfer the two-point function can be described by the 
free quark and gluon fields and thus, the integral over the hadronic spectral density equals the 
integral over the perturbative result from a certain threshold on: 



p"(s)e(s 



4) 



p 



QCD 



{s)Q{s-so) 



(2.50) 
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This ansatz will be discussed more intensively in the next section. Inserting (|2.5flj) in (|2.49|) . we 
get the final expression for the B-meson decay constant: 



4 r2 



3mi 
8^ 



so 



ds 



(m; 



2M2 



mt 



1 1 1^ \ 



2M2 J 



mliqq) 



167r 1 



+ ml 



mt 



II IT 



ml 



4M2 12M4 / 



asm) 



(2.511 



This sum rule was the starting point of many analyses through the last twenty years - see 
jl4l 1171 [3] as examples. However, the results obtained differ significantly - the estimates range 
from Jb = 130 MeV to /b = 270 MeV. The reason for this discrepancy lies mainly in the 
sensitivity of the sum rule (|2.51|) to the chosen quark mass rrib. If one takes the pole mass 
mb^poie and varies it from nih^poie = 4.6 GeV to nib^poie = 4.8 GeV, the result changes from 
Jb = 210 AleV to /b = 150 MeV. Therefore, an accurate determination of the heavy quark 
mass is necessary^. The value for the b-quark mass is usually taken from lattice calculations or 
bottonium sum rules. In [3], Jamin and Lange use the MS mass instead of the pole mass. This 
seems to improve the convergence of higher order corrections in the strong coupling a^. The 
difference of the two analyses of Khodjamirian [S] and Jamin |3] will be discussed in section IXTl 
We will give arguments in favor of choosing the MS mass and adopt the sum rule from Jamin 
and Lange as the starting point of our analysis. 



2.2 Quark-Hadron Duality 

The uncertainties in the parameters mi,, sq, M'^ and the appropriate scale, where the sum rules 
should be analyzed, are not the only limitations to the accuracy of the sum rule results. Sum 
rules have an intrinsic uncertainty. On one side, they are limited by the truncation of the op- 
erator product expansion. At higher dimensions of the expansion in l/Q^ small size instantons 
come into play and OPE cannot take these hard non-perturbative corrections into account. Fur- 
thermore, higher corrections in the strong coupling to the coefficients of the condensates and 
the perturbative calculation are scarcely known. On the other side, the hadronic spectral func- 
tion, since it is not known experimentally to all energies, is approximated by the arguments of 
quark-hadron duality. While we do not refer to the former limitations further, we will discuss 
the duality approximation more intense and give arguments for two modifications of the usual 
sum rule analysis, related to duality. 

At large negative momentum transfer, the two-point correlator is believed to be described by 
the free quark and gluon fields. We can safely neglect the power expansion on the QCD-side 
and again equate the hadronic side (|2.14|) with the perturbative side (|2.40|) at large . We get 

in the limit — > oo. We took account for the subtraction terms by using Il{Q^)/Q'^ rather 
than n((5'^) as the relevant correlator. Equation (|2.52|) is known as global duality relation. The 
integrals over the hadronic spectrum and the imaginary part of the perturbative calculation equal 

^However, the definition of the pole mass suffers from an ambiguity of order Aqcd and one cannot give a more 
precise value. This wiU be explained in section WIT] 
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in the limit of — > oo. Furthermore, the two integrands should have the same asymptotics at 
s — > oo to fulfill the equation ()2.52|1 : 

p^{s) = -ImW{s), at s^oo. (2.53) 
vr 

Setting the hadronic spectral function point-wise equal to the perturbative part is known as local 
duality. It is now argued that one can find a threshold sq, where the integral over the hadronic 
function equals the integral over the perturbative one at finite Q^, leading to the semi-local 
duality relation: 

gh so 

Thus, the uncertainty of the hadronic spectral function is replaced by the introduction of one 
more parameter. This is certainly only approximative and one has to argue for the right choice 
of the threshold-parameter sq. However, Borel transforming 1)2.54(1 introduces an exponential 
factor, which suppresses the contribution at higher momenta. After the transformation, p. 54(1 
is set in (|2.49|) and leads to the final sum rule (|2.51|) . 

In section [3.31 we will introduce a Borel mass dependent threshold sq. We can find functions 
sq{M'^) that will make fs constant over the Borel mass in a certain window. Furthermore, 
the derivative log E(M2, so(M2)), where ^(M^, so(M2)) is the right hand side of (^311) . 

can be set constant over and equal to the squared meson mass m\. The threshold sq{M'^) is 
not a unique function - we can find a continuum of functions, leading to a continuum of results 
for fs- However, we can give arguments for a preferred value of fs- 

We will also apply a second method of analyzing the sum rules, which is related to the local 
duality approach used several times by Radyushkin and collaborators [HI [7] . Taking the Borel 
parameter to infinity justifies the truncation of the condensate expansion. However at the 
same time, the exponential suppression of the contribution of the hadronic spectral function van- 
ishes. In these approaches, the threshold sq is usually taken to be around the midpoint between 
the ground state and the first resonance of the corresponding hadron. In the limit — > oo, 
the sum rule results show a rather high sensitivity on the chosen threshold sq. In section 
we suggest setting the daughter sum rule equal to the meson mass, for which a distinct value 
for So and therefore for the decay constant fs can be found. 

In the next section, we will introduce the light-cone sum rules, which will be used with the above 
modification in this thesis to get values for the weak form factor /^(O) and the coupling g^,^^. 



2.3 Light-Cone Sum Rules 

Considering processes of the type A ^ B + C , involving three hadrons or external currents, one 
needs to modify the original sum rules. A natural extension of the SVZ sum rules would be to 
sandwich three currents, interpolating hadrons by quark fields, between the physical vacuum, 

TABc{p,q) = -J d^xd%{0\T{jA{x)jB{0)jc{y)} |0), (2.55) 

and apply a similar procedure as for the two-point sum rules. However, this method, known as 
three-point sum rules, leads to difficulties in several applications. In section |2 we analyze the 
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weak form factor /^(O) for the semileptonic decay of the B-meson, at maximal recoil {q^ = 0). 
It is defined by the matrix element for the B ^ tt transition: 

{■K{p^)\ujf,d\B{pB)) = {PB +P7r)tif^{q'^) + {PB -P7r)tif'{q^) 

= 2f+{q')p., + {f^{q') + f-{q'))q,, (2-56) 

q = pb —p-K being the momentum transfer to the leptons. In the following, we concentrate only 
on the form factor f~^{q^) in front of the pion-momentum. In the case of large recoil {q^ = 0), 
the b-quark of the meson decays into a u-quark bearing a large energy of about = mb/2. 
This process is shown in figure ESI The u-quark has to recombine with the d-quark to form the 
pion. This can be done by exchanging a hard gluon over a small separation of the Fock state 
partons (hard contribution). If such an exchange is not present, the end-point regions, where 
one of the partons carries almost all momentum of the pion are enhanced (soft contribution). 

Expanding the corresponding three-point amplitude in terms of local condensates immediately 
leads to problems of convergence. Replacing the d-quark propagator by the quark-condensate 
gives contributions to the correlation function only at the end-point regions, since the conden- 
sates are made out of slowly varying vacuum fields and do not carry high momenta. This leads 
to the fact that the condensates in the OPE are accompanied by increasing powers of the heavy 
quark mass mi, |18j . One has to take into account a whole chain of condensates of increasing 
dimension to make the result numerically stable. 

This is, where the light-cone wave functions come into play |19( I2()j. Here, the expansion runs 
in the transverse distance of the partons instead of the short distance expansion of the OPE. 
This leads to a partial resummation of operators of any dimension into a single wave-function of 
given twist, where the twist of an expression is defined by the difference of its mass dimension 
and its Lorentz spin. 

Starting point for the light-cone sum rules for the weak form factor f^{(f') is the T-product of a 
pseudoscalar current, interpolating the B-meson, and a vector current sandwiched between the 
vacuum and the pion state: 

n^(9,P.) = i jd''xe''i^TT{p^)\T{u{x)j^b{x)M0)il5dm\0), (2.57) 

where q is the momentum carried away by the W-boson. To lowest order, we contract the b- 
quark fields and plug in the free propagator. This process is depicted in figure [2^ We get the 
expression 

Mq^Pn) = [m5(^(pO|n(x)7^75^^(0)|0) 

+A:"(7r(p,)|u(x)7^7,75a!(0)|0)] . (2.58) 

The two matrix elements are now expanded around = 0. The first term of the first element 
defines the twist 2 pion light-cone wave function: 

1 

(7r(p^)|n(x)7^75d(0)|0) = -ip^J^j du ^^{u)e'''P-^ + • • • , (2.59) 



where the ellipses stand for higher terms in the expansion. (pn{u) is the two-particle Fock wave- 
function of the pion - its argument u is the momentum fraction carried by one of the constituents. 
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Figure 2.6: Leading order diagram to the B 7r-transition amplitude. The u- and d-quark combine to 
form the pion state, parametrized by its distribution amphtude in the light-cone sum rule approach. 



Instead of the local condensates in the two-point sum rule approach, the correlation function 
is now expanded in non-local matrix elements, which are parametrized by the wave functions. 
Replacing in 1)2. 58() the first matrix element by (|2.59|) . we get after integration: 

1 

7 1 \2 _ — 2 ■ (2-60) 

Expanding the integrand around p^r = leads to a series in higher momenta of the wave function 
«"'(/?7r(tt), which are directly related to higher derivatives of the matrix element in (|2.59|) and 
thus to higher dimensional condensates. The crucial point is that the coefficients of this series 
are given by (at = 0): 



2n 



2 _ „2^n+l i'^2\"+l 



{ml - q 



1. (2.61) 



m 



b) 



Thus, the expansion in local condensates as one would have obtained in the 3-point sum rule 
approach would be useless, since the coefficients are not suppressed and one had to take the 
infinite series of condensates into account. In the LCSR approach this series is replaced by a 
single wave function. 

The shape of the pion wave function is not very well known. In the appendix El we give a 
definition in terms of an expansion in Gegenbauer-polynomials. This gives corrections to the 
asymptotic form (pniu) = 6u{l — u). The values of the scale-dependent coefficients are frequently 
discussed in the literature. 

The hadronic side of the light-cone sum rules for the form factor is obtained following the same 
procedure as in the SVZ sum rules for the weak decay constant. Inserting in 1)2. 57() a complete 
set of states, we get 



r \ {7T{pT:)\u-f^b\B{p^ + q)){B{p^ + q)\bijr^d\0) {TT{p^)\u-ff,b\n){n\bi-f^d\0) 



where \n) are the states of higher resonances and the continuum contribution. We rewrite these 
contribution in form of a dispersion integral over the hadronic spectral density as in (|2.14|) for 
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the decay constant. Defining the relevant correlator Il~^(q,p-^) by 

U^iq,p^) = n+(g2, (p^ + qf)p^^ + n±((?2, (p^ + qf)q^, (2.63) 
and using the definitions for the decay constant (|2.11|) and the form factors (|2.56|1 . we get 

oo 

n (■'■(>-. + ^)) = „,,(,4_(,^.^,)2) +y ,_(y, + ,)2 ^ (2-<i4) 

Applying the Borel transformation with respect to {p^^ + q)^, we finally get for the hadronic side 

n+(g^M2) = ^/^/+(g2)e-i7T + / dsp^{s)e-TI^. (2.65) 

rrib J 



Again, the hadronic spectral function will be replaced by the perturbative result via the as- 
sumption of quark-hadron duality. Thus, we need to rewrite 1)2.60(1 in form of a dispersion 
relation^ with respect to {pj^ + g)^. Taking the on-shell pion to be massless (p^ = 0), we write 
the denominator as: 

{up^ + qf-ml = u( {p^ + qf + ^J—^^ _ <\ (2.66) 



By making the substitution s = {m^ — q'^) / u + q"^ , we get 

1 oo 
n+(9 ,iq+Pn))= / , -Y = mtU / 



{s - q^){s - {p-n + qYY 
(2.67) 

Thus, we identify ipT^{s)/{s — q^) with the imaginary part of 11+ (g^, [p^^ + g)^): 

-/mn+(g2, (p^ + qf) = ^ J^^Zlif) . (2.68) 

Taking the simple duality ansatz, 

p'\s) Q{s - 4) = -/mn+(g2, (p^ + (z)^) Q{s - sq), (2.69) 
vr 

and subtracting the integral over the hadronic spectral density on both sides, we effectively cut 
the integral on the right hand side at sq. After Borel transformation with respect to [p^^ + 9)^ 
and taking back the substitution, the right hand side of the sum rule reads: 

so _ 1 22 

Khs{<l\M')=m,U j ^^^e-^=m,ul ^ ^.(n)e-''^^^^, (2.70) 

ml A 

where the lower limit in the integral is given by 

A = (2.71) 

So - r 

^Here, only a single dispersion relation is necessary, which is another advantage of the LCSR approach compared 
to the ordinary three-point sum rules. In the case of the B-meson form factor, these make use of double dispersion 
relations |21| . which leads to further uncertainties of the results. 
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Equating the hadronic and the perturbative side we get the LCSR to lowest order: 



fBr{q')e-^ = ^ / —^.{u)e . (2.72) 

Further contributions to the QCD-side come from higher terms in the expansion around the 
hght-cone = of the matrix element (|2.59|) : 

1 

{7r(pO|n(x)7^75(i(0)[0) = -ip^^U J due'^^P^'^ {^(u) + x^giiu)) 



.2 ^ 1 



+U I - ) / due''^P-^g2{u) + ■■■ . (2.73) 

^ ^0 



Here, gi and 52 ai'e two particle twist 4 wave functions. Their definition is given in the appendix 
IbI The second matrix element in (|2.58j) can be split into a symmetric and an antisymmetric 
part. The corresponding elements define the two-particle wave functions of twist 3: 

1 

{7T{p^)\u{xh5d{0m = Ufi^ jdue'^P-'^^piu), (2.74) 



multiplied by the symmetric tensor g^u. The antisymmetric matrix element is 

1 

{7T{p^)\u{x)a^ul5d{0)\0) = i{Pn^x,-p^,Xf,)^ j dne^"P-><,(n), (2.75) 



where 

= 7 ■ (2.76) 

ruu + nid 

This normalization factor is directly related to the quark condensate via the Gell-Mann-Oakes- 
Renner relation, and leads to a further error source besides the coefficients in the definitions of 
the wave functions. 



In appendixiBlwe list further contributions to the LCSR, from taking into account higher orders 
of the contraction of the b-fields. This leads effectively to an additional gluon line going from 
the b-quark into the pion wave function. The relevant amplitudes are parametrized by three 
particle twist 3 and 4 wave functions. We also list 0(as )-corrections to the twist 2 pion wave 
function. Again, all these terms have to be integrated over x and k in 1)2. 58() and then written in 
form of a dispersion relation. Here, terms proportional to higher powers of l/((up,r + 1)'^ — fri^) 
than one have to be partially integrated, which leads to further contributions to the sum rule in 
form of surface terms. 

Up to twist 4 we give the LCSR for the form factor at maximum recoil {(f' = 0) explicitly in 
chapter^ The expression is taken from Khodjamirian and Riickl [S] with first order corrections 
to the twist 2 wave function (pTr{u) taken from Bagan, Ball and Braun |22j . 
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Chapter 3 

Decay Constant 



The leptonic decay constant Jb of the B-meson is a hadronic parameter of high interest in to- 
day's physics. It appears in most calculations concerning the B-meson. Its squared value enters 
in the calculation of the B — i?-mixing amplitude, which is an important tool to measure CKM 
elements and investigate CP-violation. Therefore, a precise knowledge of Jb is strongly desired. 
However, there are no experimental measurements of the decay constant and one has to rely on 
theoretical predictions. There exist numerous sum rule analyses of the decay constant. In the 
following section, we will reanalyze the SVZ- sum rules for the decay constant to NLO-accuracy 
in the perturbative expansion, following Jamin and LangejSj. We compare the results with those 
of Khodjamirian and Ruckl[Hj. The two analyses mainly differ in the the use of the heavy quark 
mass - the running mass and the pole mass, respectively. 

We then proceed with our own analysis using the sum rule from ^ up to first order corrections 
in the strong coupling and neglecting the light quark masses. 

3.1 Discussion of Sum Rule Analyses from the Literature 

The Borel sum rules in the form of equation (|2.51j) are known for over twenty years. Unfor- 
tunately, the sum rule is very sensitive to the input parameters like the quark mass and the 
values of the threshold sq and the Borel mass M^. Therefore, the results of the various calcula- 
tions vary in a rather wide range. In 1983, Aliev and Eletskii used this sum rule to obtain 
fs = 130 MeV. This rather low value, compared to today's results, is due to the use of a quite 
high value of the pole mass (mf, = 4.8 GeV) and low value of the QCD-scale (A = 100 MeV). 
Using the same sum rule (|2.5ip . Narison extracted a typical value of about /b = 185 MeV 
with a very high threshold sq = 50 GeV"^ and a low Borel parameter ~ 2.8 GeV^ . 



In the late 90's analyses were applied using the running mass of the heavy quark in the MS- 
scheme instead of the pole mass. Their definitions are given in appendix Although the pole 
mass is independent of the renormalization scheme, it is sensitive to long distance dvnamics|23j. 
This can be observed by considering the self energy diagram. Replacing the gluon propagator 
by the full propagator with fermion loops, leads to a IR singularity of the integral. Therefore, 
an ambiguity, which is of order Kqcd^ enters the definition of the pole mass and one cannot 
assign a value to it, which is more precise. In the case of the B-meson decay constant, using the 
pole mass in ()2.51|1 . the contributions from the first order and second order corrections are of 
the same size as the one loop perturbative results, whereas the expansion is under good control, 
if the running mass is used - see [3] and also table 13. l1 
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We will now briefly demonstrate the analysis of the sum rule for fs , following Jamin and Lange 
[3] . Second order corrections to the perturbative result were found by Chetyrkin and Steinhauser 
[H] and are available as mathematica file rvs.m. Jamin and Lange also gave an expansion in the 
light quark mass up to O(m^). These corrections are numerically negligible, however, from terms 
of order O(m^) on, mass logarithms of the form log(^^/m^) appear. These can be cancelled 
by rewriting the normal ordered condensates into non-normal ordered ^31- With this at hand, 
Jamin and Lange found corrections to the quark-condensate to first order in the strong coupling 
a.,: 



"0 



ShM\so,^i) = 2 



1+1 



M2 



log ^ + - 
ra^ 3 



(3.1) 



This expression is already Borel transformed. The dependence on the threshold sq in the upper 
limit of the integral appears, when the initial expression is written as a dispersion relation and 
the continuum contribution from the hadronic side is subtracted. 

The calculation of the simple quark loop involved the replacement of the quark propagators 
with delta functions, thus setting the quarks on-shell. It follows that the heavy quark mass mf,, 
appearing in the perturbative part of the sum rule (j'i.Slj) . is per definition the pole mass. To 
express the decay constant Jb in terms of the running mass fn{iJ,), one has to apply equation 
()A.12|) . Effectively, this replacement only changes the size of the different loop-contributions - 
the total contribution should stay the same. However, the zero order perturbative contribution 
to the decay constant rises significantly after insertion of the running mass, while the first and 
second order show a convergent behavior. Thus, with the change to the running mass one can 
truncate the expansion even after the 0(as)-corrections. Keeping the pole mass leads to first- 
and second-order corrections of about the same size as the zeroth order. The truncation of the 
expansion in the strong coupling after the second order terms might lead to further errors, when 
contributions of higher order terms are not taken into account. The different contributions to 
the final results in the two schemes are shown in table (|3.1j) . 

The quark masses appearing in the leading order coefficients of the condensates are not defined 
to be the pole- or the running mass, either. We use the running mass without imposing 0{as)- 
corrections to the coefficients. 



The final sum rule for the decay constant up to 0{a1) and without corrections in the light quark 
mass reads: 
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where the right hand side is given by 
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(3.3) 



27 M2 4M2 12M4 J 

Here, A^^^ and A^^) ^re the additional first and second order corrections coming from the re- 
placement of the pole masses by use of equation (|A.12|) . In the expression (|3.3|) . the mass mi, 
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is now understood as being the scale dependent running mass fnh{p). The quark mass in the 
lower limit of the integral is kept to the pole mass. It indicates the start of the branch cut 
on the positive g^-axis. Replacing it by the running mass would lead to corrections starting at 
order 0{a^). The numerical results, if the running mass is used instead of the pole mass differ 
about 5%. A similar sum rule, except for the 0{as)- corrections to the coefficient of the quark 
condensate, was also used by Narison, who was further able to extract the running mass mfc(mb) 
at the same time |241 125| ^. 

Besides the QCD-parameters as,rnh{fi) and the phenomenological condensates, the function 
S(M^,so) also depends on the artificial parameters sq, coming from the crude duality ansatz 
H2.5U() , and the Borel mass Af ^ . Ideally, if one knew the hadronic spectral density to all energies 
and the whole perturbative expansion and non-perturbative corrections to it on the right hand 
side of the sum rules, the resulting function for the decay constant would not depend on the 
momentum transfer and thus not on the Borel mass. However, in the real world one has to 
give arguments for the values of these two parameters, where the decay constant is read off. 
After the Borel transformation, the contribution of the hadronic continuum is exponentially 
suppressed and the Borel parameter should be small enough to ensure that errors coming 
from the approximation (|2.5U|) do not contribute to a large extend. On the other side, the Borel 
mass should be high enough, such that the expansion in the condensates converges sufficiently 
fast and and the truncation of the series is justified. In our analysis of the sum rule 1)3.2(1 . we will 
take the window 4 GeV^ < < 8 GeV'^, where on the left boundary the contribution from the 
mixed quark gluon condensate contributes to less than 8% and the four quark contribution is 
negligible. The continuum contribution of the hadronic side is rather large in the whole window 
and rises to about the same size as the final sum rule on the right boundary. Therefore, the 
threshold parameter sq has to be chosen carefully. It is considered to be roughly in the range of 
the first exited state of the hadron under investigation and taken such, that the decay constant 
does not vary strongly in the Borel window. 

An often suggested procedure to get the reading point (Mq ,so) of Jb is to look at the daughter 
sum rule of (|3.3jl fI^ , defined by 

5R(m2, so) = log S(m2, so). (3.4) 

Applying the derivative on the left hand side of equation (|3.2j) it is then argued that the resulting 
function is independent of Jb and gives the squared pseudoscalar mass as a function of M^. 
Setting this function equal to the experimental mass m^, gives reading points for the threshold 
and the Borel mass. However, in our opinion this method sometimes is overrated in the literature. 
The crucial point is that the sum rule ()3.2() leads to a rather complicated function for the decay 
constant Jb, which is by no means independent of M^. The parameter that is constant in the 
sum rule (|3.2I) is the mass tub of the pseudoscalar meson. Its input value in the final sum rule 
is the experimental mass, independent of M^. Taking the derivative of the logarithm of both 
sides in (|3.2I) yields: 

ml + ^^^^^ = ^iM\so). (3.5) 

Thus, setting the ratio 5J(M^, sq) equal to the squared mass at a point Mq is just a rather 
complicated procedure to find an extremum of the initial function fB{M'^). The argument for 
taking this distinct value for the Borel mass is then not the fact that it represents the right 

^Narison adjusted the parameters sq, and the quark mass nib such that the sum rule is in the most stable 
region. He extracted mb{fnb) — (4.05 ± 0.06) GeV, which is rather small compared to typical values. 
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value for the meson mass (this is always the case) , but rather that the function for the decay con- 
stant is in a stable region around this point. In fjeure lXTl we plotted the decay constant /^(M^) 
for several values of the threshold sq and the corresponding ratios 3ft^/^(M^, sq), defined in (|3.4I) 
as functions of M^. The intersections with the experimentally measured mass itib = 5.279 GeV 
give the values of the squared Borel parameter M-^, where fs reaches an extremum. 

There is a distinct pair of the threshold and the Borel parameter, where the ratio 5?^/^(M^, sq) 
reaches a minimum at the meson mass. Taking one more derivative with respect to M"^ in (|3.5|) . 
gives 

mV|"(m2) 



+ 



fUAP) fUM^) fl{M^) 

Setting the minimum of 5R(M^, sq) equal to the mass m?^ yields the condition 

/^(M2,So) =0. 



(3.6) 



(3.7) 



Therefore, the pair (M^,so) found in this way gives a stable region, where the first two derivatives 
of the decay constant vanish, and one can extract a value for fs- However, there is no physical 
argument except the stability of the sum rule around this point that this value resembles the 
true decay constant at best. 

Using the sum rule (|3.3|) with additional corrections in the light quark mass up to order 0(mg), 
Jamin and Lange extracted for the B-meson decay constant 



Jb = (210 ± 19) MeV. 



(3.8) 



Here, the main uncertainties come from the value of the heavy quark mass rrn, = 4.21 ±0.05, the 
scale n = (3-6)Gel/ and the value of the quark condensate {qq){2 GeV) = -(267±17)3 MeV'^. 

We repeated the analysis using the procedure described above, however, only including first 
order corrections. Replacing the pole mass by the running quark mass leads to corrections to 
the 0(as)-terms of: 



Ad) 



mf 



3 , 

2 + -log^ 
2 



(3.9) 



For the running b-quark mass, we take the value from El-Khadra j27j . mf, = 4.24 ± 0.11 GeV, 
which is obtained as an average from a large collection of different analyses, including lattice 
simulations and sum rule methods for several systems containing b-quarks. The corresponding 
pole mass is nih^poie = 4.69 GeV, using the two-loop relation (|A.12|) . The values for the other 
input parameters we used in this analysis and the following ones are listed in table l3?2l Taking 
the mean values we obtained 



Jb = 208 MeV 



so = 34.1 GeV^ 



= 5.5 GeV^. 



(3.10) 
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43468 


33064 


1855 
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-4203 


-449 



f2 

J B,KR 
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/!" 


.(0) 


/(-) 


32314 


12045 


11170 


11388 


-2290 



Table 3.1: Contributions of the different terms to the squared decay constant fs for the two analyses of 
Jamin and Lange (left table) and Khodjamirian and Riickl (right table). indicates the contribution 
from condensates with mass dimension four or higher. 
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Figure 3.1: Decay constant (graph 
a) and the Ratio SR^/^ (b), defined 
in equation H3.5|) . as functions of the 
Borel parameter Af^, plotted for three 
different values of the threshold sq. 
The intersections of K^^'^ with the 
experimentally measured meson mass 
niB = 5.279 GeV correspond to ex- 
tremal points of the decay constant. In 
the case of sq = 34.1 GeV, also the 
second derivative of the corresponding 
M'^[GeV^] function of the decay constant vanishes. 







b) 




so = 34.9 GeV3^ 








SQ = 34.1 GeV'^^^^^ 


s„ = 33.7 GeV2~\„____^ ^^^^.^^^"^ 
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Although we only included first order corrections and used a different quark running mass and 
corresponding pole mass, this result is very close to ()3.8|) . We do not give an estimate of the 
error, since we are only interested in the relative size of the different contributions to the decay 
constant. 

The contributions to from the several terms of the sum rule are shown in table 13.11 This 
table shows also our result for the sum rule analysis following Khodjamirian and Riickl By 
using the suggested values of mi,^poie = 4.7 GeV for the b-quark pole mass and sq = 35 GeV"^ for 
the threshold, we could exactly reproduce the value Jb = 180 MeV . Note that the first order 
corrections in the strong coupling are of about the same size than the zeroth order perturbative 
result and contribute to about 35% of the final result. By switching to the running mass the 
main contribution of the perturbative calculation is shifted to the zeroth-order term, whereas 
the 0(as)-corrections are about 4% of the final result. Including the second order corrections 
from Chetyrkin and Steinhauser jH] in the pole mass scheme, the contribution is again of the 
same size than the zeroth and first order corrections. In the analysis of Jamin and Lange the 
0(a^ )-corrections are negligible. 

Since the expansion in the strong coupling constant is under good control compared to the sum 
rule using the pole mass for the heavy quark, we will use the sum rule in the form of equation 
(|3.2|1 up to 0{as) as the starting point of our own analysis. 

3.2 The Limit ^ oo 

As argued above, the method of setting the daughter sum rule ()3.4() equal to the meson mass is 
equivalent to search for a threshold parameter sq, for which the sum rule for the decay constant 
is in a stable region. Another method was suggested several times by Radyushkin [HI [7j in 
sum rule analyses of the pion form factor. By taking the limit oo, one can get rid off 
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higher dimensional condensates and the sum rule is again in a stable region. In this limit, the 
uncertainties coming from the truncation of the power series vanish, whereas higher resonances 
and the continuum contribution are not suppressed exponentially anymore. Thus, the duality 
approximation is the crucial assumption in this approach and the source of the main error. In 
the limit — > oo, the sum rule reads: 



4 f2 



8^ 



so 



ds 



{rn\ 



b.poLe 



+ml 



TT 



- + log — 

b 



nib 



+ h-GG) 

II TT 



(3.11) 



Apart from the constant contributions of the two condensates, the decay constant is now de- 
scribed as the average of the transition of a free quark pair of invariant mass s in the duality 
interval rnl^^i^ < s < sq. In the analysis of the pion form factor the threshold is estimated 
to be the midpoint between the first (vr) and second (Ai) resonance. Applied to the B-meson 
sum rule, the threshold should then be around sq ~ 31 GeV"^, assuming that the mass of the 
first excited state B' is about 0.5 GeV higher than the ground state. 



Instead of setting the threshold to the above value and varying it within a certain interval, we 
propose to take a distinct value of sq. In the limit — > oo, the sum rule for the decay constant 
is only a function of the threshold. Its first derivative scales like oc Although Jb is 

asymptotic for any value of the threshold, we suggest taking the value Sq, that satisfies 

mVI'(m2,.S)*'^°°0, (3.12) 

for which the sum rule reaches its asymptotic value the fastest. Using (|3.12|) in equation (|3.5|) . Sg 
can be found by setting the daughter sum rule of 1)3. 3|) to the meson mass in the limit — > oo: 

?Sl{M^,s*Q)^'^°^ml. (3.13) 

In figure OISl we demonstrate the matching of 5R(M^, sq) to the meson mass for different values 
of the running quark mass and show the functions of the corresponding sum rules for the decay 
constant /b- Using the values of table we find a mean threshold of Sq = 31.7 GeV. The 
corresponding decay constant is: 

/b = (178 ± 18 ± 20) MeF. (3.14) 

Here, the first error is inflicted by the uncertainties in the quark mass. We separate this error 
from the remaining ones, since when the sum rule for the form factor in section [4.11 is divided 



Parameter 


Value 


AfB 


Taken from 




(4.24 ±0.11) GeV 


±18 




im) 


-(267 ± 16)3 MeV 


±5 


m 




(3.5 - 6) GeV 


±14 




i^GG) 


(0.012 ±0.012) GeV^ 


±1 


IHI 


ml 


(0.8 ±0.2) GeV"^ 


±0 


u 




(8 ±8) -5 GeV^ 


±0 





Table 3.2: Values of the parameters used for the analysis oo and the corresponding errors they 

inflict on the final value for the decay constant fs — 178 MeV. 
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mj = 4.13 GeV 



mt = 4.24 GeV 
mt = 4.35 GeV 



""0 10 ZO 30 JO bO 10 ZO 30 HO bO 

Figure 3.2: a: Matching of tlie daugliter sum rule (|3.4(l to tlie meson mass iriB — 5.279 GeV in the Hmit 
— > oo for different values of the running quark mass mf,. - b: The corresponding sum rules for the 
decay constant, using the threshold Sq, found by the matching procedure in the left figure. 




by fs and the quark mass is varied simultaneously, the uncertainties coming from the mass will 
almost cancel in the final result for f'^(0). 

The uncertainties of the parameters used and the corresponding errors to the final result are 
listed in table l3?2l The decay constant is very sensitive to the running quark mass and the scale. 
At lower values of //, the perturbative expansion becomes quite important. The inclusion of sec- 
ond order corrections in the strong coupling would give a better control of the scale dependence. 
However, we do not intend to give a too accurate result, for the accuracy of this method itself is 
limited intrinsically. Since the exponential damping of the higher excitations and the continuum 
contribution is taken away by the limiting procedure — > oo, errors coming from the crude 
duality ansatz are rather emphasized. We therefore sum up linearly all errors inflicted by the 
several input parameters except the quark mass, leading to the relatively large overall error in 
(Tmi) . 

Although the decay constant fs is very sensitive to the values of the input parameters, the 
corresponding threshold Sq is not. Varying the parameters within the errors indicated in table 
13.21 and adjusting the threshold such that the daughter sum rule matches the meson mass 
in the limit oo (|3.13j) . Sq always stays in the rather small interval 

31.5 GeV"^ < So < 32.0 GeV'^. (3.15) 

Vice versa, the value of the decay constant is quite sensitive to the threshold sq. Thus, taking 
the limit — > oo and guessing the threshold to be in a certain region, instead of fixing it by 
the method described above, leads to large uncertainties of the final result. 

In section 13.41 we will relate the value 1)3. 14() with the result of the analysis of the next section 
and compare our results with other analyses from the literature. 

3.3 Borel Mass dependent Threshold so{M^) 

In this section, we propose a further approach to analyze the sum rule (|2.51|) . As argued above 
and illustrated in figure l3Tl the method of setting the logarithmic derivative (|3.4j) equal to the 
meson mass is a rather sophisticated way to search for a value of the threshold sq, for that the 
sum rule of the decay constant is in a stable region. This works nicely for the B-meson decay 
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constant - the resulting function of /b is almost constant over the Borel window. However, there 
are sum rules, as will be seen later, in which the resulting function for the hadronic parameter 
has a strong dependence on M^. Here, one cannot find a threshold sq in a physical accessible 
region for which the first two derivatives of the sum rule vanish and sometimes not even the first 
one. 

Ideally, the decay constant is not a function of the Borel parameter. From a certain value of 
on, the expansion in the vacuum condensates in equation (|.S.3j) is under good control and 
its truncation does not induce a significant error to the final result. Assuming that all errors to 
the sum rule come from the duality approximation (|2.54)) . we introduce a Borel mass dependent 
threshold so(M^) that compensates this error. This means, we search for functions so(M^), for 
which the decay constant fs is not a function of the Borel mass. 

At first sight, this modification is quite similar to the initial approach of adjusting the constant 
threshold such, that the sum rule is in stable region. Instead of having a decay constant that 
is dependent on M^, this dependence is now shifted to the threshold ■so(Af^). Furthermore one 
does not find a unique value of the decay constant. For each value of Jb, one should be able to 
find a function so(-/Vf^) . However, this approach has an intuitive interpretation. It is rather clear 
that the hadronic spectral density can only be represented approximately by the perturbative 
free quark spectral function. Introducing in (|2.54|) a g^-dependent threshold can compensate for 
this error. For each distinct value of the momentum transfer, one can than impose a value s{q'^) 
rendering the approximation exact. After Borel transformation, the resulting sum rule for the 
decay constant should then be independent of provided that the errors induced from the 
truncation of the power series and the expansion in the strong coupling are negligible^. 

It still remains the question of where to read off the value of the decay constant, since one can 
find a continuum of functions so{M'^), for which the sum rule is constant over M^. At this 
point, the duality assumption reenters the analysis. Assuming (|2.54|) is a good approximation 
in the sense that one can find a constant threshold for which the two sides equal, we suggest 
looking for the most stable function so(M^) to reflect this duality ansatz. Thus, whereas in most 
other applications the threshold is a constant and its value enters by physical reasoning or is 
found by searching for the most stable function /(M^), the proposed modification features con- 
stant functions /b in a given Borel window and we search for the most stable functions so(M^). 

Since the sum rule for /^(M^), which results from setting the minimum of the daughter sum 
rule equal to the meson mass is already quite constant in the window 4 GeV < < 8 GeV 
(see figure m|) . it is not surprising, that the following analysis will lead to the same result. 

Taking the derivative M^d/dM'^ of the logarithm of both sides of equation (|3.2j) yields: 



Setting the right hand side equal to the squared meson mass m^, we get a differential equation 
for the function so(M^): 



■^In the analysis below, we introduce the functional dependence of the threshold after the Borel transformation, 
since the limiting procedure of equation 12.471 would be non-trivial if one included the dependence on the momen- 
tum transfer before the transformation. Nevertheless, the physical picture, that shifting the g^(M'^)-dependence 
to the threshold is compensating errors coming from the duality approximation, stays the same 



m% + 



M^fl'jM^) ^ M4gf^S(M^go(M^)) + M^4(M^)^S(M^^o(M^)) 
/2(M2) - S(M2,so(M2)) 



(3.16) 




(3.17) 
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Figure 3.3: Sum rule analysis for Borel mass dependent thresholds. Graph a) shows various functions 
for so{M'^), which are solutions to the differential equation (|3.17() . Below, the corresponding functions 
for the decay constant are plotted (b). Graph c) shows the variation of the threshold functions so{AP) 
plotted in a), as a function of their average value in the given Borel window. Diagram d) shows the decay 
constant plotted over the average values of the thresholds. 



We do not write the explicit form of the right hand side, since it is a rather complicated function 
of so(M^) and its derivation is straightforward. Setting the solutions of equation H3.17|) into the 
sum rule for the decay constant will yield constant functions of fs- 

This equation was solved numerically. We started at certain values for the threshold sq at 
= 6 GeV^ and the functions were then evaluated with mathematica, using a numerical 
procedure with adaptive step size. 

Figure 13. 31 a shows some solutions of so(M^) we obtained for the mean values of the input 
parameters listed in table The corresponding decay constants are given below (b). As 
argued above, we read off the decay constant at the most stable function of so(M^). This could 
be done by simply looking at the graphs of the threshold so(-/Vf^), since the method itself contains 
implicit errors and one should not seek for a too accurate result. However, to demonstrate the 
equivalence to the result (|3.1()|) obtained by following the analysis of Jamin and Lange [HI, we 
introduce a variation of the threshold functions: 

8 

s:ar = \JisUM')-^fdM'. (3.18) 
4 

In figure EiSlc, we listed the variations of the functions so{M'^) over their average value. Graph 
I3.3l d shows the corresponding decay constants, again plotted as a function of the mean value of 
the threshold. Reading off the decay constant at the point with the least variation gives 

fB = 208 MeV. (3.19) 
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Figure 3.4: Borel mass dependent functions for tlie tlircshold so(Af^) and corresponding decay constants 
/b, plotted over the average value of the thresholds. On the left side rat, ~ 4.35 GeV was used as input 
value for the running quark mass - on the right side rrifc = 4.13 GeV . 



Thus, we reproduced the result (|3.10|) from the analysis in section I3.H following Jamin and 
Lange. 

In table 13.31 we listed the value of the input parameter we used for this analysis and the 
corresponding errors they inflict on the final result. In figure im we plotted the results of our 
analysis for the b-quark masses rrih = 4.13 GeV and mf, = 4.35 GeV . Again, Jb is very sensitive 
to the running mass and the scale. We again separate the error coming from the quark mass 
from the remaining errors, added in quadrature: 

fB = (208^23 MeV (3.20) 

The total error of fs is rather large. This is mainly due to the large uncertainties, which we 
assigned to the quark mass. In the following two chapters, we will analyze the sum rules for the 
B-meson semi-leptonic form factor /^(O) and the coupling g^'B-K- These sum rules are propor- 
tional to the decay constant. When the two sum rule results are divided by /b and the quark 



Parameter 


Value 


-Mb 


+A/b 


Taken from 




(4.24 ±0.11) GeV 


231 


167 


123 


im) 


-(267 ± 16)3 MeV 


200 


212 


m 




(3.5 - 6) GeV 


192 


221 




i^GG) 


(0.012 ± 0.012) GeV^ 


207 


208 


i 


ml 


(0.8 ± 0.2) GeV^ 


208 


203 






(8±8)-5 GeV'^ 


208 


208 



Table 3.3: Input parameters used for the analysis with a Borel mass dependent threshold. The column 
— A/b(+A/b) shows the deviations from the central result Jb = 208 MeV , when the lower(higher) limit 
of the input parameter is used. 
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masses are varied simultaneously, the final errors due to the uncertainties in the quark mass 
decrease drastically. 

We also checked that changes in the boundaries of the Borel window have a minor effect on 
the corresponding variations of the functions so{M'^) and the point of extraction of the decay 
constant. 

3.4 Summary of the Results and Relation to the Literature 

The two results pT^ and differ about 15%. In the sum rules we analyze in this thesis, 

the method of taking the limit — > oo will always lead to significantly lower results than 
the method of imposing a Borel mass dependent threshold so(M^). Considering only the LO- 
perturbative term of the daughter sum rule 

so 

J dssf{s)e~^, (3.21) 

ml 

it is clear from the positivity of the integrand f(s) and the vanishing of the exponential damping 
in the limit oo, that the relevant value of the integral needed in the fitting procedure 

is reached at a rather small threshold sq. In the second method, we are looking for functions 
So(M^) that do not depend strongly on the Borel mass. Thus, the contribution from the term 
proportional to Sq(M'^) in (|3.16|) is small compared to 1)3.21(1 and therefore, the upper limit 
so(M^) of the integral has to be higher than in the former case. The difference in the size of the 
two thresholds translates into the discrepancy of the final results for the decay constant. This 
can be inferred from the figures and IXTl Since the limit 3?(M^,Sq) m?^ is reached from 
below, the corresponding sum rule for the decay constant is monotonously decreasing and in the 
Borel window of figure ITT] it would lie below the plotted curves. 

In table ITU we listed our results for the decay constant and several results from the literature. 
The latest sum rule analyses jH] and j25j . including three-loop corrections to the perturbative 
calculations and using the running quark mass, coincide with our result from using a Borel mass 
dependent threshold. The result from Wittig |2iij is a global estimate of recent unquenched lat- 
tice calculations. It is about 10% lower than the sum rule results. Nevertheless, all the quoted 
values are in agreement, when the error bars are taken into account. 



fB[MeV] 


Ref. 


Method 


178 ± 27 
2081^11 
180 ± 30 

190 ±9 
210 ± 19 
205 ± 23 
173 ± 13+f 
218 ± btli 

191 ± 23^?9 


This Thesis 
This Thesis 
Khodjamirian '98 8 
Narison '98(24j 

Jamin '010 
Narison '01^3] 
Abada APE'99[2H1 
Bowler UKQCD'00|2ni 
Wittig '03[1]1 


SVZ-SR 0{as) 
SVZ-SR 0{as) 
SVZ-SR 0{as) 
SVZ-SR 0{a,) 
SVZ-SR 0{al) 
SVZ-SR 0{al) 
Lattice 
Lattice 

Lattice-summa 


running mass; lim Af ^ — > cxd 
running mass; so(M^) 
pole mass 
running mass 
running mass 
running mass 

ry of unquenched results 



Table 3.4: Comparison of our results with a collection of values from other analyses in the literature. 
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Before we proceed with the analysis of the B- meson weak form factor /^(O) and the strong 
couphng gB*B-Ki present briefly the analysis of the two point sum rule for the vector meson 
decay constant fs* ■ 



3.5 Weak Decay Constant / 



B* 



The decay constant fs* is needed in the analysis for the strong coupling Qb^b-k-: since the sum 
rule we use to extract this coupling is proportional to the decay constants of the B- and the 
S*-meson. We will analyze the light cone sum rule for the coupling by the use of the two new 
methods introduced in this chapter. To get a consistent value for Qbob-k^ ^'^^ divide the sum 
rule by the corresponding decay constants Jb and fs* , obtained within the same method. We 
therefore summarize in this section our analysis of the B* decay constant. 

We collected the sum rule for the heavy- light vector meson constant from Dominguez j31j and 
Reinders ,12]. Neglecting the numerically insignificant contribution of the four quark condensate, 
we get up to mass dimension 5 in the condensates and first order in the perturbative expansion: 
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where the first order correction in the strong coupling is given by: 
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(3.22) 



bs"^ — m^s — 2m^ 
(s - ml){2s + ml) ' 



(3.23) 



We again replaced the pole mass in the integrand by the running MS'-mass. This lead to the 
additional term: 

6ml{s + ml) f^^3^^^_fj^ 



A«(^ 



4 I 2 

m^ + 



2s2 



mt 



(3.24) 



As before, with the pole mass replaced by the running mass, the sum rule shows better conver- 
gence. The first order correction amounts to about 20% of the perturbative part, whereas before 
this contribution was about 60%. The use of the running mass increases the overall perturbative 
contribution about 30%. 



Parameter 


Value 


-Mb- 


+AfB* 


Taken from 


mb{mb) 


(4.24 ±0.11) GeV 


221 


174 


123 




-(267 ± 16)3 MeV 


192 


203 






(3.5 - 5) GeV 


182 


208 






(0.012 ±0.012) GeV^ 


197 


198 


i 


ml 


(0.8 ± 0.2) GeV^ 


197 


197 






Table 3.5: Values of the input parameters used in the sum rule analysis of fs* in the limit AP oo 
and the corresponding errors they inflict on the final result. The central result is /b' = 197 MeV. 
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Parameter 


Value 




±A/b* 


Taken from 


mhirrih) 


(4.24 ±0.11) GeV 


281 


213 






-(267 ± 16)3 MeV 


233 


257 


m 




(3.5 - 5) GeV 


226 


258 


- 




(0.012 ±0.012) GeV^ 


245 


246 




mo 


(0.8 ±0.2) GeV"^ 


252 


243 





Table 3.6: Values of the input parameters used in the sum rule analysis of fs* with a Borel mass 
dependent threshold sqIM"^) and the corresponding errors they inflict on the final result. The central 
result of this analysis is /s* = 245 MeV. 

We analyzed this sum rule using the two methods introduced in the preceding sections. The 
values for the decay constant we extracted are: 

/b. = (197 ± 24 ± 16) MeV, (3.25) 

when the limit ^ co is taken and 

fs* = (245 ± 35 ± 23) MeV (3.26) 

in the case of a Borel mass dependent threshold so(M^). Again, the first errors come from 
the uncertain running quark mass rnb{fnb) = (4.24 ± 0.11) GeV. The second error summarizes 
further uncertainties. The individual contributions of the uncertainties of the input parameters 
and the dependence on the scale ^ are collected in the tables 11^31 and 

The result (IT^ is close to a NLO-result from Khodjamirian|321, fs* = (195 ± 35) MeV, 
whereas ()3.26() is far off that value. This can be traced back to the strong dependence of the 
sum rule on the Borel parameter if a constant threshold sq is used. To stabilize the sum rule 
by introduction of a Borel mass dependent threshold will require rather high average values of 
so(M^). Nevertheless, we will use the two decay constants from this analysis in the extraction 
of the couphng ^b^b^- 
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Chapter 4 

B-Form Factor /+(0) 



In this chapter, we analyze the weak semileptonic form factor /^(O) of the B-meson at maximal 
recoil. In section l2?3l we already gave the sum rule for the form factor of the twist two pion wave 
function. Up to twist 4 we adopt the sum rule from Khodjamirian |S]. For = this is: 



/B/+(0)e- 



4 



2m2 



B 



(4.1) 



where ^{so,M'^) is given by: 



$(so,M2) 



du ""b 
u 



fwiu) + — I U(pp{u) + 

mi, \ 6 



1 + 



mt 



2u]VP 



Ami 



Ami 



9i{u) + 



dtg2{t) 1 + 



mt 



uA'P J 



(4.2) 



Here, the lower limit in the integral is now given by A = m^/so- A~^(so, M'^) is the surface 
term, coming from the partial integration necessary to build up a dispersion relation. The ex- 
plicit form of the wave functions as well as the three particle contributions and the first order 
corrections to the twist 2 wave function is given in appendix IbI 



The right hand side of equation (|4.2)) is effectively an expansion in powers of 1/ ttM^ rather than 
just M^. Taking the Borel window of the sum rule for the decay constant /b, 4 GeV^ < < 
8 GeV"^, and dividing the boundaries by an average value of the integration variable u, we obtain 
a new window that takes the change from the parameter to uM^ into account. Evaluating 
the integral in (|4.2|) in the Borel window for the decay constant, we find an average value of u 
varying between 0.5 and 0.7, for which the integrand is equal to the evaluated integral. Thus, 
we use for the analysis of the form factor a window with boundaries of about the doubled values 
as for the decay constant: 

8 GeV^ < A'P < 14 GeV"^. (4.3) 

We checked that the hadronic contributions of higher states and the continuum, which are sub- 
tracted from the QCD side by the duality assumption are of about 25% of the final result on 
the right edge of the Borel window. Higher twist contributions are under good control and 
the contributions of the three particle wave functions are almost negligible. For central input 
parameters, we depicted the contributions from the single terms to the total sum rule in figure 
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Figure 4.1: Contributions to the final sum rule for /b/^(0): twist 2 (full line); 2P - twist 3 (dashed); 
0(as)-twist 2 (dash-dotted); 3P - twist 3 and twist 4 (dotted). We used the central values mf, = 4.7 GeV 
and So = 35 GeV^ as input parameters. The sum of all contributions is plotted in figure ^21 



14.11 The dotted line collects the two particle twist 4 and all the three particle twist 3 and 
4 contributions. The sum of these terms is only about 2%. Therefore, the truncation of the 
expansion after these terms is justified. 



We will evaluate the sum rule at the scale 

/xfo = yn4^^ = 2.4 GeV. (4.4) 

This scale, which is often used in analyses of hadronic properties of the B-meson, reflects the 
average virtuality of the b-quark in the meson. A further ambiguity lies in the use of the quark 
mass. Nevertheless, taking the running mass mb(mb) = 4.24 it 0.11 GeV and running it down 
to the scale (|4.4|) . we get mb{fib) = 4.69 GeV, which is of about the same size as the pole mass 
obtained from the two-loop relation HA.12I) . Therefore, we will take the typically used value of 
the pole mass: 

mb = A.7±0.1GeV (4.5) 

The errors are also of the same size as for the running mass. In figure IX2k . we plotted /b/^(0) 
for the three different masses rrib = {4.6 GeV, 4.7 GeV, 4.8 GeV} and sq = 35 GeV^. Other in- 
put parameters and the values of the coefficients of the wave functions can be found in appendix 

m 



Taking the mean values = 11 GeV^ and sq = 35 GeV^, we read off /s/+(0) = 0.049 GeV, 
consistent with the analysis of Khodjamirian and Riickl, who extracted f~^{0) = 0.27 with a 
rather low central value for the decay constant, /b = 180 MeV. However, the sum rule is quite 
unstable in the given Borel window as can be seen in figure It decreases monotonously in 
about 10%. In fi gure 14. 21 b. we plotted the corresponding daughter sum rules. They are 
clearly lower than the experimentally measured meson mass niB = 5.279 GeV. Intersections 
with the mass of the B-meson will only occur at higher values of the Borel parameter or the 
threshold sq. A unique pair (M^,so) that would correspond to a vanishing second derivative 
of fBf^{0){M'^) is never reached within finite values of the Borel mass and the threshold. The 
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Figure 4.2: a: Sum rule for the product /b/^(0) for sq = 35 GeV and quark masses mi, — 
(4.6, 4.7, 4.8) GeV. Graph b) shows the corresponding daughter sum rule. The curves do not intersect 
with the meson mass, and therefore /s/^(0) decreases monotonously and does not reach a minimum in 
the Borel window. 

choice of the read off pair (M^, sq) is therefore rather ambiguous. This motivates the apphcation 
of the two modified analyses introduced before. 

4.1 The Limit ^ oo 

In this section we will apply the method of local duality, introduced in section 13.21 to the sum 
rule for /b/~^(0). Taking the limit — > oo justifies the truncation of the power expansion in 
and minimizes uncertainties coming from higher twist wave functions. On the other side, 
uncertainties in the hadronic spectral function are emphasized, since they are not exponentially 
damped anymore. 

We again take the derivative d/dM'^ of the logarithm of the right hand side in equation (|4.2|) . 
We then adjust the threshold sq such, that the square root of this daughter sum rule takes on 
the value of the meson mass uib = 5.279 GeV in the limit — > oo. Using the threshold 
obtained in this way in the sum rule 1)4. the product fsf^iM"^) reaches its asymptotic value 
very fast. In figure UHl we plotted the matching of the daughter sum rules for the center value 
of the quark mass and its largest errors m^y = (4.7 it 0.1) GeV , as well as the corresponding 
functions for fsf^i^) of M^. The central value is 

lim fBf^{M^) = 0.0463 GeV, (4.6) 

and the error inflicted by the uncertainties of the quark mass is about 12%. 

Since the decay constant obtained with this method was rather sensitive to the scale /x, at which 
the sum rule was evaluated, we also calculated the values of fsf^iO) at the scales /i = 2 GeV 
and /i = 3 GeV. There are many //-dependent parameters in 1)4. 2() and we will give their scaling 
laws in the following. 

The scaling behavior of /Xtt can be obtained from the Gell-Mann-Oakes-Renner relation (|2.42j) and 
is therefore directly related to the running of quark masses. The scaling law of the coefficients of 
the twist two pion wave function iPtt{u) is given in ()B.2|) . The scale dependence of the parameters 
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Figure 4.3: Matching of tlie daugiiter sum rule to tlie experimental meson mass ttib — 5.279 GeV in 
the limit — *■ oo (a). The values for the thresholds obtained by the matching are used to determine 
the values for the product /b/^(0) in the same limit (b). 



appearing in the twist 3 and 4 wave functions are taken from Braun and Filianovj33j: 



a 



/stt 



"s(/ii)y 

Iflzt 

2 



a 



/37r'^l,0l/^lj 



a. 



2 / 122 



Stt 



5 

5 45 

for the parameters of the two and three particle wave functions of twist 3 and 



4a;2,o 



(4.7) 
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a 



(4.10) 
(4.11) 



for the twist 4 wave functions. Using these scahng laws, we evaluated the sum rule (|4.1j) at 
the scales /i = 2 GeV and /x = 3 GeV . The two functions of the daughter sum rule are almost 
similar and the matching leads to two thresholds, which are also very close to each other: 



so(2 GeV) = 35.46 GeV, so(3 GeV) = 35.45 GeV. 



(4.12) 



Setting these two thresholds in the corresponding functions for fsf^iM"^), we get a small scale 
dependence of about ±4%, when varying the scale fi = 2 GeV to /i = 3 GeV. This is shown in 
figure |^3]a. We finally obtain: 



lim /b/+(m2 



(0.046 ± 0.006 ± 0.002) GeV, 



(4.13) 



where the first error is due to uncertainties in the quark mass and the second comes from the 
scale at which the sum rule 1)4. 2(1 is evaluated. There are further errors coming from the un- 
certainties of the input parameters fi-j^ and /sjt, however, they are small compared to the values 
above and we do not take them into account in this analysis. 
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Figure 4.4: a: Sum rule for /b/+(0) evaluated at the scales = 2 GeV and fi — 3 GeV. b: /b/^(0), 
when the asymptotic pion wave function is used (dashed line), compared to the central result (solid line). 



A rather large error is inflicted by the uncertainties of the coefficients of the light cone wave 
functions, which are not known to a high accuracy. We therefore also analyzed the sum rule 
for the asymptotic twist 2 pion wave function. In figure U^Jj, we plotted the sum rule for the 
asymptotic wave function, compared to our central result. The value for the product fsf^i^) 
decreases about 13%. This number can be interpreted as a rather conservative error estimate 
on the result (|4.13j) . 

For consistency, we divide (|4.13() by the decay constant (|3.14() . obtained with the same method 
in section 13.21 In the ratio, the two errors from the quark masses almost cancel, when they 
are varied simultaneously. Thus, although both sum rules show a high sensitivity to the chosen 
quark mass, this dependence becomes very weak if the two results are combined. We obtain: 

/+(0) = 0.26 ±0.05. (4.14) 

The error due to the unknown quark mass amounts to less than 4%. We included the error 
estimate of 13% from the uncertainties of the coefficients of the wave functions. Added to the 
other error sources, the final error estimate is about 20%. In section ESI we relate our results 
to values found in the literature. 

The two individual analyses of the decay constant Jb in section l3?2l and the product fsf^i^) 
in this section show rather low values compared to other methods. However, since the ratio of 
the two values is taken, some errors cancel as shown above and the error of about 20% of the 
final result is reasonable. It is hoped that the intrinsic error coming from the method itself also 
shrinks when the ratio of the two results is taken. This will be highly supported in the next 
section. 



4.2 Borel Mass dependent Threshold so(M^) 

The sum rule for the form factor (|4.1|) leads to functions that have a rather strong dependence 
on for typical values of the input parameters, as can be seen from figure U3 With higher 
values of the thresholds the functions develop a minimum in the Borel window, however, the 
second derivative does not vanish as was the case in the analysis of the decay constant. Thus, 
one cannot find a unique pair (M^,so) for which the daughter sum rule is equal to the meson 
mass at its extremal point. 
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Figure 4.5: Sum rule analysis for fsf^iO) with a Borel mass dependent threshold so{]VP). Diagram 
a) shows several solutions of the differential equation H4.15|l . Plot d) shows the corresponding values 
of /b/^(0) over the mean values of the thresholds. Graph c) shows the variations of the functions of 
the threshold to give the region where the value of /s/^(0) should be read off. Diagram b) shows 
the functions /b/^(0) when the solutions so(Af^) are inserted and gives a crosscheck of the numerical 
analysis. 



In this section we apply the method outlined in section I3..SI to the sum rule of the form factor. 
Thus, we shift the functional dependence of the form factor to the threshold sq by solving the 
differential equation 

'° M4^$(M2,.o(M2)) 

where ^{M'^ , so{M^)) is given in 1)4. 2|) . For central values of the input parameters, we plot- 
ted several solutions in figure 14.51 We again plotted the variations of the curves to read off 
the value for the product /b/^(0) at the most stable threshold. This threshold is rather high 
(so ~ 40.4 GeV"^) compared to the earlier analyses. As can be seen in figure the dependence 
of the sum rule on the Borel parameter is rather strong for smaller values of the threshold. Thus, 
one needs to impose a function sq{M'^) that is also strongly varying with the Borel mass. To 
higher thresholds, this dependence of the sum rules on becomes weaker, which is then also 
the case for so{M'^). 

In figure EH we plotted the results from using mf, = 4.6 GeV and = 4.8 GeV as the values 
of the quark mass. As before, the error to the product fsf^iO) coming from the variation of 
the b-quark mass is large. It is about 15% of the final result. 

We also examined the scale dependence of /b/^(0) extracted from this analysis. Like before, 
we evaluated the sum rule at /i = 2 GeV and fj, = 3 GeV. The results are shown in figure ETTI 
The values of fsf^iO) are quite stable to variations of the scale /i. The error is less than 4%. 
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Figure 4.6: Analyses for the different values of the quark mass: mb = 4.8 Gey(left side) and rrib = 
4.6 GeV{ngh.t side). The upper diagrams show the solutions for the thresholds so(M^) - the lower 
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Figure 4.7: Scale dependence of the sum rule for fsf'^iO)- On the left side, the results are shown for 
the scale iJ = 2 GeV and on the right side, the scale fj, = 3 GeV was used. 
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Figure 4.8: Sum rule analysis for the asymptotic pion wave function. The product /s/^(0) decreases 
significantly about 15%. 

We again analyzed the sum rule using the asymptotic pion wave function to take the error from 
the uncertainties of the coefficients of the wave functions into account. The resulting thresholds 
so{M'^) and values for fsf^ifi) can be seen in figure 

Neglecting further error sources, we extract for the product /b/'''(0): 

fBf^iO) = (0.055 ± 0.009 ± 0.008) GeV, (4.16) 

where the first error is the one coming from the variations of the quark mass and the second is 
due to the scale dependence and uncertainties of the coefficients of the wave functions. If we 
divide this result by the decay constant 1)3. 20() . that was analyzed using the same method, we 
obtain: 

/+(0) = 0.26 ±0.05. (4.17) 

The change of the quark mass to the extremal values in the two sum rule results (|4.16|) and 
p.20() leads to slightly higher ratios (about 3%) than the central value in either case. Again, the 
main error comes from the conservative error estimate of the uncertainties in the coefficients of 
the wave functions. 

4.3 Summary of the Results and Relation to the Literature 

It is amazing that the two results 1)4. 14() and 1)4. 17() are so close to each other. The central results 
differ by less than 2%. The two methods have very different physical meanings. Whereas the 
first method of taking the limit — > oo emphasizes errors coming from the duality approxima- 
tion, the second method tries to suppress these uncertainties. Therefore, the discrepancy of the 
individual results for the decay constant and the product /s/"^ was rather natural. It was hoped 
that when the ratio of the two sum rule results was taken, the intrinsic errors from the methods 
would be lowered. However, we did not expect that the two results for the semi-leptonic form 
factor are that close to each other. In the next chapter we will analyze the sum rule for the 
strong coupling 9^*^^- The sum rule is proportional to fs and Jb* and will be divided by these 
constants, analyzed within the same method. However, this time, the extracted couplings differ 
about 13%. 

We did not take 0(as )-corrections to the two particle twist 3 wave functions into account. These 
were recently obtained in ^35j. The 0(as )-corrections to the twist 2 wave function raised this 
contribution about 30%. Assuming that the new corrections are also about 30% of the twist 3 
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/+(o) 


Ref. 


Method 


0.26 ± 0.05 
0.26 ± 0.05 
0.24 ±0.025 
0.27 ±0.05 
0.28 ±0.05 

0.26 ± 0.06 ± 0.05 
0.28 ± 0.04 


This Thesis 
This Thesis 
BaU '91 HU 
Khodjamirian '97|S1 
Bagan '97|^ 

BaU '01 Eg 
Abada APE'99|2B1 


LCSR; 0(as)-twist 2 corrections; hmM^ oo 
LCSR; 0(as)-twist 2 corrections; so{M'^) 
3-pt SR 

LCSR; 0(as)-twist 2 corrections; 

LCSR; 0(as)-twist 2 corrections; average of full and 

asymptotic twist 2 pion wave function 

LCSR 0(Q!s)-twist 2 and 3 corrections 

Lattice 



Table 4.1: Comparison of our results to a collection of values from the literature 



wave functions, we expect the final result of /^(O) to raise about 10%. 

In table lO we collected some values for the form factor f~^{0) from the literature. Our results 
seem to fit quite well in this table, especially when the large error bars are taken into account. 
Furthermore, the results from LCSR analyses and 3-point sum rules coincide with lattice cal- 
culations. This will be very different in case of the strong coupling g^*^^ analyzed in the next 
section. 
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Chapter 5 

Coupling g^^^ 



In this chapter, we will apply our two methods on the sum rule for the strong coupling constant 
9b'B-k- '^^^ process B* Bn is kinematically forbidden, thus one has to rely on theoretical 
predictions. This coupling is closely related to the B ^ ir form factor. It is assumed that the 
form factor is dominated by the vector meson at zero recoil. In the single pole model, the form 
factor at large ~ m^, is parametrized by 

Thus, the coupling is part of the residue of the pole of the form factor. It was observed in the 
literature that the vector meson dominance is also valid to smaller values of the transferred 
momentum. However, at = 0, the two approaches, the pole model and the LCSR for the form 
factor, differ in their scaling behavior in the heavy quark limit rrib — > oo (see |36j). We therefore 
do not compare our results for the form factor at maximum recoil with the sum rule results for 
the coupling gB*Bn- 

In the next section, we will briefly sketch the derivation of the sum rule, following 36^ and then 
apply the modified analyses to it. 

5.1 Derivation of the Sum Rule 

We take the following definition for the i3*i37r-coupling: 

{B*\q)7r'ip^)\B-{p^ + q)) = -^^.^^eX' (5-2) 

where is the polarization vector of the i?*-meson. Matrix elements of states with different 
charges are related to the definition of the coupling by isospin symmetry. 

Starting point is again the correlation function (|2.57jl of two pseudoscalar currents, sandwiched 
between a pion state and the vacuum. This function can be split in two parts 

n^(p^, q) = n(g2, (p^ + qf)p^^ + n(g^ (p^ + qf)q^, (5.3) 

and again, we are only interested in the function in front of the pion momentum. 

On the hadronic side, we first insert a set of intermediate states carrying B-meson quantum 
numbers. This leads to the product of two matrix elements: 

{TT\uj^b\B){B\bij5d\0) (5.4) 
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The second element is proportional to the B-meson decay constant. Using crossing symmetry 
and inserting another set of intermediate states, this time carrying i?*-meson quantum numbers, 
we obtain the matrix element of the upper definition H5.2() and a second element, which defines 
the vector-meson decay constant: 



(5.5) 



Using this definition, we analyzed the sum rule for fs* ■ The results were summarized in section 
13:51 



The hadronic side of the correlation function is: 

_ m^niB'fBfB* 93*37, 



+ 



mb{q'^ - m%,){{p^ + q)"^ - m%) J {si - q^){s2 - {pn + q)"^) 



P^{si,S2)dsidS2 



+ 



(5.6) 



The ellipses stand for subtraction terms, which will be removed when the Borel transformation 
is applied. The limit S/j of the double dispersion integral is now a region in the (si, S2)-plane. 

When the expressions (|5.6)) and the corresponding QCD-side are Borel transformed with respect 
to q'^ and (pt^ + g)^, the two sides are equaled and the hadronic continuum is subtracted by the 
assumption of quark hadron duality: 



4 



So 



fBfB*9B*Bne 1 e 



rrib 



^QCD^ ITfg ^^dsids2. 



(5.7) 



The derivation of the QCD spectral function is exactly the same as in section [2121 since we used 
the same initial correlation function. The explicit subtraction procedure is rather complicated 
and we refer to the appendix in 36 . The two Borel parameters are believed to be of the same 
magnitude. The choice = Ml = 2M^ simplifies the derivation of the final expression for the 
sum rule: 



fBfB*gB*BTT(^ 



ml 



-r(M2,so), 



(5.8) 



with r(M^,so) given by: 



r(M2,so) 



Mm + ^ f ^v^.(i/2) + ^^.(1/2)) + 

rrib \2 3 / m;,/^ 



+e m2 



V9,(l/2) +52(1/2) 



4m^ 
M2" 



1/2 

5i(l/2) + j g2{v)dv^ +4p^ 




2sn 



+ 



37r 



/ 2 
mf 



2 e dt. 



(5.9) 



2ml 



Here, the three particle twist 3 and 4 contributions are given by: 



i: 



(3) 



i: 



3P 
(4) 



3P 



32 



(20 + 5uJifi + wi,! - 4u;2,( 



(5.10) 
(5.11) 
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Figure 5.1: Sum rule for /s/^.g^.^^ and corresponding daughter sum rule. The threshold is set to 
So = 35 GeV and the graphs are plotted for rat, — (4.6,4.7,4.8) GeV . The sum rule decreases about 10% 
in the Borel window. 



The 0(as)-corrections to the two particle twist 2 wave function were obtained by Khodjamirian 
and collaborators, and we adopt the expression from [32] : 



, vr" , /x\ , / x\ 3(3j;3 + 22x^ + 40x + 24) , /x 
m = -+31og(-)log(l + -)-J ilog(- 



2(2 + x 

- 3Li2{-x) - 3Li2{-x - ' 
3(3j;2 + 20x + 20) 6x{l + x) log(l + x 



+6Li2 i^-j - 3Li2i-x) - 3Li2{-x - 1) - 31og(l + x) log(2 + x) 



4(2 + x)2 + (2 + x)3 ■ ^^'^^^ 

For the following analyses, we took the same window for the Borel parameter as for the form 
factor in chapter^ 

8 GeV < < 14 GeV^. (5.13) 

The contribution of the twist 2 wave function at the left boundary is about 60% and of the 
twist 3 terms about 40% of the total. The twist 4 contribution is numerically negligible. The 
contribution from the hadronic continuum grows from 20% to about 30% at the right edge of 
the window. Thus, the use of this window is reasonable. 

For typical values of the input parameters - sq = 35 GeV, rub = (4.7 ± 0.1) GeV - we plotted in 
figure the sum rule for the product fBfB*gB*B-K corresponding daughter sum rule, 

defined by: 

^{M\so) = M^^logT{M\so). (5.14) 

Applying the same derivative on the left hand side of equation ()5.8() . we see that the sum rule 
reaches a minimum, when the daughter sum rule is equal to 

= ^ 28.11 GeV^ (5.15) 

if the experimental values for the meson masses tub = 5.279 GeV and tub* = 5.325 GeV are 
taken. The daughter sum rule is rather far apart from rh and the sum rule for the product 
fBfB*gB*B-K quite unstable in the chosen window. 

In the following two sections, we will apply the methods introduced in the preceding two chapters 
on the sum rule (j5.8j) . 
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Figure 5.2: Matching procedure for the sura rule of fEjE'dB'S-n- Diagram a) shows the matching of the 
daughter sum rule in the limit A/^ oo for three different values of the quark mass mf,. The values of 
the threshold obtained by this procedure are used to extract the value for /b/b* g^.^.^. in same limit (b). 



5.2 The Limit ^ oo 

In this section, we will analyze the sum rule 1)5. 8(1 by imposing local duality. We take the limit 
oo and match the daughter sum rule (|5.14j) to the mass m, defined in equation (|5.15|) . 
From this matching, we can extract a unique value for the threshold sq and therefore a unique 
value for fBfB*9B*B7v 

In figureEH we plotted the graphs of the sum rule and its daughter sum rule for the central value 
of the quark mass nib = 4.7 GeV and the two extremal values = 4.6 GeV and mb = 4.8 GeV. 
The error due to the uncertainties in the quark mass are about 7%. Besides this error source, 
we also analyzed the scale dependence of the sum rule. In figure l5^ a. we plotted the resulting 
graphs for ^ = 2 GeV and /x = 3 GeV. Again, we took account for the uncertainties of the 
coefficients of the wave functions by analyzing the sum rule for the asymptotic twist 2 wave 
function (figure Einb)- The resulting deviation from the central result (9%) was assigned as 
additional error to the final result. Neglecting further error sources, we extract: 

fBfB*9B'B. = (0.71 ± 0.05 ± 0.07) GeV^, (5.16) 

where the first error comes from the uncertain quark mass and the second summarizes the error 
from the scale dependence and the uncertainties of the coefficients of the wave functions. The 
value (|5.16|) is about 10% higher than the result obtained in (HHI, which did not include 0[as)- 
corrections, and again about 10% smaller than the result from |32j . which was the first analysis 
using the 0(as)-corrections to the twist 2 wave function. 

If we divide 1)5. 16() by the corresponding value of the decay constant 1)3. 14() . extracted using the 
same method, we can get the residue of the form factor in the simple pole model: 

c = ^^*^B*B^ = (0.37 ± 0.06) GeV. (5.17) 

zniB* 

Here, we simultaneously varied the quark mass, leading to a rather small contribution to the 
overall error (about 4%). This value is compatible with the result from the NLO-analysis |32j - 
(0.42 ± 0.09) GeV. 

Dividing (|5.16j) by the pseudoscalar decay constant (|3.14|) and the vector-meson decay constant 
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Figure 5.3: a: Scale dependence of the sum rule for the product fBfB'gB'Bir ™ ^^'^ limit A'P oo. 
When the scale is varied from ji = 2 GeV to = 3 GeV , the resulting value of the sum rule decreases 
about 7%. - b: Comparison of the sum rule result using the asymptotic pion wave function (dashed line) 
with the central result. The coupling Ob'Btt increases about 9%, when the asymptotic function is used. 



p.25|) . we get for the vr-coupling: 

= 20±5. (5.18) 

This time, the uncertainty of the quark mass inflicts an error of about the same size as the 
further error sources to the final result, which is about 17%. The result (|5.18|) is noticeably 
lower than the result of the LO-analysis ^36 {9^*^11 = 29 it 3) and of about the same size as the 
NLO-result ^ - c/^*^^ = 22 ± 7. 

In section l5^ we will discuss the results for the coupling gB*BiT and relate them to the literature. 

5.3 Borel Mass dependent Threshold so(M^) 

In this section we analyze the sum rule (|5.8j) with the method introduced in section 13.31 We 
take the threshold as a function of the Borel mass so(M^) and solve the differential equation 

rn^^(M^ .o(M^)) - M^^^(M^ so{M^)) 

to get constant functions for the product fBfB*gB*B-K- Since the initial sum rule was rather 
dependent on the Borel mass (see figure l5?T1) . we expect the resulting functions so(M^) to have 
quite high average values. 

For central values of the input parameters, the solutions so(M^) and the corresponding sum 
rules are shown in figure 15.41 The average values of the threshold are quite high compared to 
the values in the last section. This is equivalent to the behavior of the initial sum rule, which 
gets more stable with higher values of the constant threshold sq. 

In figure Em we plotted the functions ■so(M^) for the two extremal input values for the quark 
mass rrib = 4.6 GeV and = 4.8 GeV . Whereas we could find a minimum of the variation 
of the threshold functions in the case of = 4.8 GeV , we could not do so for m^y = 4.6 GeV . 
Solving the differential equation (|5.19j) leads to functions so(M^), which are strongly dependent 
on the Borel mass (see figure 15. 51 c) . Their variations are rather high compared to the former 
analyses and do not show a minimum. Nevertheless, we could find a region where the variations 
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Figure 5.4: Sum rule analysis of the product fsfB'gB'BTT- Graph a) shows solutions to the differential 
equation H5.19|l . b) the corresponding functions for fBfs'gB'Bn' shows these values plotted as functions 
of the average threshold so- Graph c) shows the variations of the thresholds s(M^). The product 
/B/B*<?B»B7r is extracted at the minimum. 



are changing very slowly. We plotted the corresponding values fBfB*gB*Bn ^^^^ region in 
figure l^3b . However, to be on the safe side we added the difference of the central value and the 
lower bound from rrib = 4.8 GeV to the mean result to obtain an upper bound. 

The scale dependence of this sum rule analysis is shown in figure ESI and is less than 4% of the 
final result, if the scales /i = 2 GeV and /i = 3 GeV are used to extract fBfB*9B*B-K- 

We also analyzed the sum rule for the asymptotic twist 2 wave function. The results are shown 
in figure 1^771 As in the last section, the central result increases about 9%. 

The final result for the product fBfB*gB*B-K extract, is: 

fBfB*9B'B. = (0.88 ± 0.08 ± 0.09) GeV^ . (5.20) 

This is about 10% higher than the result obtained by Khodjamirian and collaborators using the 
same sum rule P^ . 

Using this result to again extract a value for the residue c in the single pole model, we get after 
dividing by 

c= (0.40+^) Gey, (5.21) 
which is close to the result 1)5. 17() obtained in the preceding section. 

Dividing (|5.20|) by the two decay constants (|.S.20|) and (|3.26j) . we get: 

Qb^b. = 17^^- (5.22) 
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Figure 5.5: Sum rule analyses for different values of the quark mass. The left side shows the thresholds 
So(M^) and corresponding values plotted over so(M^) for the case nib = 4.8 GeV. The right 

side summarize the analysis for nib = 4.6 GeV. 
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Figure 5.7: Sum rule analysis for the asymptotic twist 2 pion wave function. The extracted value is 
about 9% higher than the central result. 

The difference between tlie positive and negative error assigned to tlie results 1)5.22(1 and (|5.21|) 
is due to the different effect of a higher and lower quark mass on the decay constant (|.S.20|) . The 
error inflicted on the final results by the different masses is about 18%. 

5.4 Discussion and Relation to the Literature 

In the last two sections, we analyzed the sum rule for the product /b/b'^^'Btt- chapter 
1^ it was hoped that when this sum rule is divided by the corresponding sum rule results for the 
decay constants, the intrinsic errors of the two methods would be reduced. The results (|5.17|) 
and ((5.2111 for the residue c of the pole in the simple pole model ((5.1(1 are within 6% of each 
other, which is a quite satisfying, considering the significant difference of the two initial sum 
rule results ((5.16(1 and ((5.20(1 . However, the extracted values for the coupling 5^.^^ ((5.18(1 and 
((5.22(1 differ about 15%. The reason for this deviation might be the unusual large value for the 
vector meson decay constant, found in section ESI 

We did not take further perturbative corrections into account. Assuming that the first order 
corrections to the twist 3 contribution is about the same relative size as the 0(as )-twist 2 cor- 
rections (~ 40%), we get another uncertainty of 20%. 

Furthermore, the hadronic continuum contribution was only subtracted from the leading terms 
in the expansion ((5.9() . The inflicted error to the final result is negligible as well as variations 



Sb*Btt 


Ref. 


Method 


20 ± 5 


This Thesis 


LCSR; 0(as)-twist 2; limM^ ^ oo 


ntl 


This Thesis 


LCSR; 0{as)-twist 2; so(M^) 


15 ±4 


Colangelo '94|S3 


NLO-soft pion limit 


28 ±6 


Belyaev '95^ 


LCSR-LO 


14 ±4 


Dosch '96|SHI 


Double Moment SR 


22 ± 7 


Khodjamirian '99(5^ 


LCSR; 0(a^)-twist 2; 


14.5 ±3.9 


Navarra '00 39^ 


3-pt SR 


32 ±5 


Melikhov '0102 


Dispersion Approach 


47 ± 5 ± 8 


Abada 'OSEH 


Lattice, quenched 



Table 5.1: Comparison of our results to a collection from the literature. 
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from the boundary Sq in the plane (si, S2) |32] . 



In table I^TTl we listed our results for the coupling and a collection from the literature. It should 
be noted that the values obtained in sum rule approaches are significantly lower than results 
from quark models and lattice calculations. In the case of the D-meson coupling g^)*^^^, the 
situation is quite similar. Here, a first measurement of the coupling gave a value almost twice as 
large as the sum rule results [l^]. Several ideas were proposed to cure this problem of the sum 
rules |43| I44j. We will discuss this issue in the conclusion of this thesis. 
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Chapter 6 

Conclusion and Outlook 



In this thesis we introduced two modifications of the analysis of SVZ sum rules and light-cone 
sum rules. These were aimed to extract the relevant hadronic property in a stable region of the 
sum rule. However, the two methods differ significantly in their physical meaning. In the first 
method, we took the limit oo. In this limit, the exponential suppression of the higher 

states and the continuum vanishes. Therefore, this method rather emphasizes uncertainties from 
the duality ansatz. The result becomes quite sensitive to the chosen threshold sq. We fixed this 
threshold by setting the daughter sum rule equal to the meson mass in the limit — > oo. In 
the second method, we imposed a Borel mass dependence of the threshold and we chose it such, 
that the sum rules became constant in a given Borel window. With this procedure, we tried to 
damp the uncertainties of the duality approximation. The sum rule result was extracted from 
a stable region of the threshold functions so(M^), arguing that the duality assumption is still a 
good approximation and the corresponding thresholds So(M^) should not vary strongly. 

The individual sum rule results of the two methods gave rather different values of the hadronic 
properties. The values extracted from the limiting procedure oo were always significantly 

lower than the values from the second method. However, when the ratio of two sum rule results 
was taken, as was the case for the B-meson form factor and the strong coupling, the intrinsic 
errors seemed to compensate partially. The final results of the form factor were almost identical. 
The results for the coupling came closer to each other, compared to the initial sum rule results 
for the product f 3/3*9 b-'Btv 

In table El we listed the results from all analyses in this thesis. We also gave the average value 
for each hadronic property. Since one of the methods is emphasizing errors from the duality 
approximation and the other is trying to suppress them, it seems reasonable to give the mean 
value of the results obtained with the two methods, since the true value might lie in between 
these two extremal points of view. 

We already related our results to values from the literature. Whereas the extracted values for 
the decay constants and the form factor seem to fit quite well to other results, the value for 
the coupling §3*3^ is close to other sum rule results and far off from values obtained within 
other models and lattice calculations. In the case of the D-meson coupling, an experimental 
measurement confirmed the values from quark models and lattice calculations. However, the 
experimental value was about 40% higher than the sum rule results^. 

'^The value l|5.20|l extracted from imposing a Borel mass dependent threshold so(M^) is about 20% higher than 
the value one would extract from figure IF!T1 and the inclusion of 0(as )-corrections to the twist 3 wave functions 
could raise the result by another 20%. However, when the result H5.2()|l is divided by the decay constants extracted 
from the same sum rule analysis, we get a rather low value for the coupling. 
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lim oo 




average result 


fB[MeV] 
fB*[MeV] 
/+(0) 
9b*b-k 


178 ± 27 
197 ± 29 
0.26 ±0.05 
20 ± 5 


208+^^ 
245 ± 42 
26 ±0.05 

1713 


193+^^ 
221 ± 26 
26 ±4 

19^3 



Table 6.1: Summary of the hadronic properties extracted from the two different sum rule methods. In 
the last column, we gave the corresponding value when the average of the two sum rule results is taken. 



One of the main uncertainties of sum rules is the scarcely known hadronic spectral function. 
Nevertheless, the crude ansatz of taking the first resonance and relating the remaining contin- 
uum to the QCD-side of the sum rules gave many results in agreement with experiment. In the 
case of the strong coupling the initial sum rule is very unstable in the Borel window (see figure 
15.11) . In 133]) Becirevic et al. suggested to take the first radial excitations of the D- and the 
Z)*-meson explicitly into the calculations. This introduces several further parameters. However, 
they were able to show that the inclusion of a negative radial excitation improves the stability 
of the sum rule. Furthermore, the region of best stability of the sum rule coincides with the 
experimentally measured value. It would therefore be interesting to incorporate the radial exci- 
tation in the analyses of this thesis. Ideally, one would expect a significant raise in the extracted 
coupling g^*^^ and minor changes to the decay constants and the form factor. 

In this thesis, we analyzed hadronic properties of the B-meson. These analyses can easily be 
applied to the D-meson. However, evaluating the sum rules at scales of about = 1 ~ 2 GeV 
demands a very careful treatment of perturbative corrections. The inclusion of higher order 
corrections should improve the stability with respect to variations of the scale /j,. In the light of 
the coming experimental data at CLEO we try to give predictions in the near future. 
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Appendix A 

Renormalization Group Properties 



Any Green's function in QCD has to obey the CaUan-Symanzik equation to make it independent 
of the artificial scale-parameter that comes in the calculations during the renormalization pro- 
cedure. The solution of the Callan-Symanzik equation implies a scale-dependence of the strong 
coupling constant as{p). The runnning of the coupling is described by the Gell-Mann-Low 
/3-function, which can be written as a power series of the coupling: 



_ das 



diogK 

Mo 

The leading coefficients of this series are: 



= -6n-^ - h 
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bo 
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(A.2) 



The third coefficient is scheme dependent and we have given it in the MS'-scheme. n/ is the 
number of quark flavors below the scale /x. Solving for the strong coupling up to two loop order 
yields: 



as(/x) = 



l + ^6olog^ 



1 - 



bias{fJ,o) 



Mo 



4Trbo + as{no)bl log 



log(l + ^^^6olog£) I (A.3) 



This equation can be parametrized by introduction of a mass scale A 



An 



bo log ^ 

To one-loop order, the scale A is defined by: 



1 



2bi log log ^ 
^0 log^ 



log 



[4 

A2 



1. 



(A.4) 



(A.5) 



The scale A depends on the number of flavors used in the calculations. In the numerical analyses 
we will use the experimental value of the strong coupling at the Z-Boson mass, taken from the 
Particle Data Group: 



as{Mz) = 0.1172(20); Mz = 91.1876(21) GeV 



(A.6) 



54 



We will then run the coupling down to the scale at which we are evaluating the sum rules. At 
the threshold of the b-quark mass, one has to apply matching conditions [461 147j . These are 
obtained by constructing an effective theory with = 4 light quarks and one heavy quark and 
matching it with the full theory at the b-quark threshold. The matching is done by considering 
several Green's functions and require consistency up to a certain order in (l/m;,) at the thresh- 
olds. Taking only continuous matching conditions of the strong coupling results in a rather 
strong dependence on the actual scale of the matching. In these cases, one usually accounts 
for the uncertainties by introducing errors via matching of the two effective couplings at 2mh 
and mh/2. We are very grateful to M. Jamin for making available a sophisticated mathematica 
package for running of the strong coupling and the quark masses. It implements the matching 
conditions from Chetyrkin|37] up to four-loop accuracy. 

Treating the quark masses like coupling constants in perturbative QCD, we get from the solution 
of the Callan-Symanzik equation a prescription for the running of the masses: 

Mo 

The 7's are the anomalous dimensions of the scalar gg-operator. They are given by 

70 = 2 

101 5 , 

71 = nt. (A. 8) 

' 12 18 ^ ^ ' 



To one-loop accuracy we get 



m\^^) = m\^^l) ( 1 + ^^ho log , (A. 



47r 



which can be rewritten as: 



'log fly' ' 







The running quark-mass is defined during the renormalization of the quark mass. It is introduced 
to cancel the divergent part of the propagator. Up to first order in the running coupling it is in 
the M5-scheme|lH]: 

rn{^i) = mo (l + ^) , (A.ll) 

where niQ is the bare quark mass and oj is the deviation from four dimensions, D = 4 — 2uj. Being 
defined in this way, the running mass is different from the pole mass, the location of the pole of 
the full renormalized quark propagator [Wj. Unlike in the leptonic case, where the pole mass of 
the fermions can be identified with their physical mass, quarks are never observed on-shell due 
to confinement. Defined in this way, the two quark masses can be related to each other: 




rripoie = rn{ii') 1 + - + log ^ + O(a^) . (A.12) 



Higher loop contributions can be found in j5()| , however the analyses in this thesis will not extend 
to this accuracy. 
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Appendix B 

Light-Cone Wave Functions 



In this appendix we give the definitions of the two- and three-particle hght-cone wave functions 
up to twist four, and hst terms, contributing to the form factor (|4.2|1 in section 0] 



B.l Definitions of the Light-Cone Wave Functions 

The twist 2 pion wave function can be expressed as a sum of orthogonal polynomials, by using 
partial wave decomposition. The derivation can be found in j5T] : 



V9^(n, n) = 6u{l -u)(l + a2{fi)C^^^{2u - 1) + a4(^)cf ^(2^ -!) + •••). (B.l) 

Cn^^ are Gegenbauer polynomials. The coefficients are not known very well. Efforts were 
made to extract them from several sources, like two-point sum rules, lattice and instanton 
physics. They are multiplicative renormalizable: 

a„(^) = aM (^4^) , (B.2) 



with 

3(n + l)(n + 2) 3 f-f 
In our analysis, we will take the values from |34j : 

a2{fib) = 0.218, a^ifih) = 0.084, (B.4) 

at the scale ni, = 2.4 GeV. 

Further definitions of two- and three-particle wave functions are taken from Khodjamirian [S]. 
The relevant references where they are collected from, can also be found in this review. For con- 
venience, we list them in the following, where we did not write explicitly the scale dependence 
of the wave functions. 

The three particle twist 3 wave function (p^.,^ and the corresponding two particle twist 3 wave 
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functions ipp and ip^, which follow from (^3,^ by equations of motion are given by: 

fSTvicii) = 180010203 (2 + ujifi{7a3 - 3) + 2u;2,o(2 - 4Q;ia2 - 803 + 803) 

+2a;i,i(3aia2-2a3 + 3ai)) (B.5) 

(^p(u) = l + i^(3(2u- 1)2-1) 

+t|^(4cj2,o - ^1,1 - 2wi,o)(35(2n - 1)^ - 30{2u - 1)^ + 3) (B.6) 

3/37r I ^ 1 . \ /c/r,,. 1 \2 



V^{u) = 6u{l-u){l + ^^i^5--oJi^o]{5{2u-lY-l) 



^(4cu2,o - c^i,i)(21(2^x - 1)^ - 14(2n - 1)^ + l) . (B.7) 



16/i„ 

Here, fij^ = 2.02 GeV is given in (|2.76|) . The other parameters at the scale ^ib are given by: 

wi,o = -2.18, 102,0 = 8.12, wi,! = -2.59, fs^ = 0.0026 GeF^. (B.8) 

The twist 4 three particle wave functions and the corresponding two-particle wave functions are 
defined as: 

99x(a,) = 3052(ai - 02)0^ Q + 2e(l- 203)^ (B.9) 

ipwicKi) = 1206'^ e{ai — a2)aia2C(3 (B.IO) 

^x(a*) = 30^2(1-03)0^ Q + 2e(l- 203)) (B.ll) 

^1 



(^11 (a,) = -12052^10203 (^3 +£(1-303) j (B.12) 

g^(u) = -5'^n-ufu^ + -e5^(u(l-u)(2 + l3u(l-u)) + 10u^logu(2-3u + -u^) 

2 2 \ 5 

+10(1 - uf log(l - u)(3n - 1 + ^(1 - uf)^ (B.13) 

92{n) = Y^\(l-u){2u-l). (B.14) 
And at the scale /Ub, we have: 

52 = 0.17 Gey2, £ = 0.36 (B.15) 

B.2 Contributions to the semi-leptonic Form Factor 

The surface term in (|4.2() is given by (at g2 = 0): 

il\ 4 so\ fml\ 4 ^51(5 



So c?"" 



ml 
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The three-particle contributions are: 



; du J dtti (5 ( 1 — 
2/3. 



a,; 



(ai + MQ3) 



1 



<^37r(ai) 1 



(ai + 'ua3)M2 



(B.17) 



While these two contributions were taken from Khodjamirian 1^, we will use the first order 
corrections to the twist 2 wave function from [^. These are easier to handle numerically and 
one can also see explicitly the soft and hard contributions to the heavy-light decay in the heavy 
quark limit rrih 00: 



/+^'^(0)e^ = -^r«(5o,M2) 



where 

t(^)(so,m2) 
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(B.18) 
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and 



Ei{x) 

Li 
L2 



If. , (t - mlf 
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(B.20) 

(B.21) 
(B.22) 
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